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general theory of wall interference 

FOR STATIC STABILITY TESTS IN CLOSED RECTANGULAR 
TEST SECTIONS AND IN GROUND EFFECT 

By Harry H. Heyson 
Langley Research Center 

SUMMARY 

A theory is developed which predicts the interference velocities and interference 
velocity gradients caused by the walls of the tunnel. Large wake deflections are allowed 
in both the lateral and vertical directions. The theory includes V/STOL and conventional 
wall-interference theories and ground effect as special cases. Symmetry and interchange 
relationships between the interference factors are developed, and extensive numerical 
results are presented. 

Use of the interference factors to correct data depends upon the availability of 
detailed aerodynamic treatment in nonuniform flow of the model under test. In most 
tests the available aerodynamic treatments will be found either inadequate or too time 
consuming for rigorous routine correction of data relating to lateral-directional stability. 

INTRODUCTION 

Reference 1 presents a review (published in 1966) of the status of subsonic wall- 
interference theory. In addition, reference 2 presents a similar review of wall- 
interference theory as it pertains to models with highly deflected wakes. In general, 
these reviews find that existing theory is reasonably adequate when used to determine 
the effect of the walls upon the overall performance of the model. Corrections to longi- 
tudinal stability measurements, obtained by calculating the interference at the tail of the 
model and by examining the gradients of interference over the wing chord, are somewhat 
less satisfactorily verified. However, provided that the model is reasonably small, 
generally satisfactory results can be obtained for longitudinal stability tests as well, pro- 
vided that considerable care is exercised not only in the wall-interference calculations 
but also in correcting for nonuniformities in the distribution of flow angularities within 
the basic wind-tunnel flow. 

Theoretical treatment of wall interference as applied to lateral-directional stability 
testing is essentially nonexistent. In a few cases, the asymmetries generated by yawing 



a finite wing have been noted by reference 1 to be significant; however, the effect of the 
side forces generated by the asymmetric model are not considered in the available 
treatments. 

A complete wall-interference theory including lateral-directional stability tests 
must consider a number of features. First, since the model is asymmetrically disposed in 
the tunnel, side forces will be present. These side forces will be influenced by the walls 
and will result in interference velocities just as do the lift and drag forces. Secondly, 
since lateral-directional testing primarily constitutes a study of the moments rather than 
the forces on the model, and since it is primarily the gradients of the interference veloci- 
ties which affect moments, it is necessary to examine all of the possible interference gra- 
dients in the tunnel. Observe that the mutual consideration of the above two features is a 
major complication, since, in order to consider all of the possible wall influences, it 
becomes necessary to examine the interference velocities in the orthogonal directions as 
caused individually by forces in these three directions for a total of nine interference 
velocities, as well as the gradients of these nine interference velocities along each of the 
three coordinates for a total of 27 interference velocity gradients. Finally, since modern 
aircraft developments have resulted in the need to test many types capable of extraordi- 
narily large wake deflections, it is desirable to examine the effect of these wake deflec- 
tions, both vertically and laterally, on the resultant interference factors. That such 
effects may be large has already been demonstrated by the more limited analyses pre- 
sented in references 3 to 8. 

The present paper presents an analysis which considers all of the foregoing features. 
In the basic theoretical treatment, the model is assumed to be vanishingly small and 
located at an arbitrary point within the tunnel test section, and the interference velocities 
and their gradients are obtained at an equally arbitrary point within the test section. The 
wake, which actually follows a curved path from the model to infinity downstream in the 
tunnel, is linearized to consist of a series of straight-line segments which follow approxi- 
mately the same path. Except for modifications to allow for side forces and for a lateral 
wake deflection, this wake is essentially the same as that used in references 5 and 6. 

The assumption that the model is vanishingly small is less restrictive than appears 
on the surface. Linear superposition of the results may be used to obtain the appropriate 
interferences for models of arbitrary size or configuration just as in references 6 and 7. 

A few sample calculations of this nature are included in the present paper. 

In the application of the present theory it is necessary to estimate the wake deflec- 
tion angles with respect to the tunnel axes. The momentum analysis presented originally 
in reference 9 is not adequate in the present case because that paper did not consider the 
possibility of side forces with respect to the coordinate system. The necessary modifica- 
tions to the theory of reference 9 are derived herein. Furthermore, references 10 and 11 
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have already noted that it is necessary to modify the wake angles as computed by momen- 
tum theory in order to account for wake rollup. A few remarks on useful approximations 
to the rolled up wake deflection angles are also included. 

As presented herein, the theory pertains specifically only to completely closed rec- 
tangular wind tunnels. Formal extension to completely open or to closed-on-bottom-only 
test sections, as in reference 5, is simple. However, it should be observed that the deri- 
vation of the usual boundary condition imposed at an open boundary depends upon the use 
of small perturbation assumptions (ref. 12). If the wake actually impinges upon an open 
boundary, these assumptions may be severely violated (refs. 3 and 5). Under such condi- 
tions the theoretically obtained interference velocities may be grossly in error. It is 
recommended that such extension should only be attempted with great caution both in the 
theoretical treatment and in the application of the results to the correction of wind-tunnel 
data. 

Similar concern must be observed when attempting to correct data from closed test 
sections if the wake deflections are sufficiently great. It was first shown experimentally 
(refs. 13 and 14) and later theoretically (refs. 8 and 11) that sufficiently large wake deflec- 
tions can result in such enormous alterations of the flow within the test section that the 
measured data no longer represent any free-air flight condition (although under certain 
circumstances they may approximate flight in ground effect (ref. 8)). Reasonably satis- 
factory correlations of the conditions under which these effects limit testing (ref. 2) have 
been obtained for tests which involve essentially no lateral wake deflection. The equiv- 
alent limiting conditions for wakes with large horizontal as well as vertical deflections 
are completely unexplored. 

The effect of wall interference on the measured data will differ for different models 
according to the sensitivity of the model characteristics to particular interference veloci- 
ties or velocity gradients. The corrections to data can be no better than the investigator's 
ability to calculate the effect of these velocities and velocity gradients on the model char- 
acteristics. For some classes of models, particularly many V/STOL types, there is an 
inadequate theoretical background with which to calculate corrections. For other types 
of models, theoretical means for these calculations may exist but be too lengthy for prac- 
tical application to large masses of data. An exploration of actual correction formulas 
would involve almost all known aerodynamics and is obviously beyond the scope of any one 
paper or any one author. No such complete treatment is intended herein. Instead, a few 
remarks on the treatment of wall effects as a problem in similitude are included in the 
hope that this discussion will point out sources of information presently available for 
approximate data corrections. 
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SYMBOLS 


In order to arrive at a set of consistent axes and signs in the present analysis, it 
has been necessary to define certain quantities in a manner in conflict with many stability 
analyses. The reader is cautioned to consider carefully the definitions, particularly of 
positive directions, and provide for suitable conversion to his desired standard. 

(f f) Sin X H Sin X V + ( ? h) Sin X V cos X H + f h) sin X H cos X V 

(? |) sin X H sin Xy - ^ sin x y cos x H + (? sin x H cos x v 

Aq reference area in ground effect, 4h^ 

A m momentum area of lifting system 

A t wind-tunnel cross-sectional area, 4BH 

ai functions relating length along the wake to the generalized coordinates xj 

(see eqs. (22) to (25)) 

B semiwidth of wind tunnel 

b lateral distance from origin of doublet wake to right-hand (viewed from 

behind) sidewall of wind tunnel 

C D drag coefficient, ^ 

Cjjy = - COS^Xjj cos^Xy 

C L lift coefficient, ^ 

C, rolling-moment coefficient, positive when moment tends to roll aircraft 

to left R o llin g moment 
’ qSc 

C m pitching-moment coefficient, positive when moment tends to pitch aircraft 

Pitching moment 
nose up, B 

qSc 

C n yawing-moment coefficient, positive when moment tends to yaw aircraft 

noS e left, Jawing moment 
qSc 


A = h J r 0 C hv - 
A' = h|R 0 C H y - 
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"R 


'M 


r> 

resultant- force coefficient, -=£ 

qS 

lateral-force coefficient, positive along Y-axis, ^ 


qS 


jet-momentum coefficient, Momentum 

qS 


mean aerodynamic chord, or equivalent dimension for nondimensionalizing 
moments 


C V 

D 

H 

h 

L 


cos 


*H 


COS Xy 

drag, force directed along X-axis 

induced drag, induced force directed along X-axis 

semiheight of wind tunnel 

height of model above floor of wind tunnel (or ground) 
lift 


distance from origin along wake 


m,n,p,q,r integers 


m< 


doublet strength for doublets with axes directed along the i-axis 


n 


ratio of final to initial induced velocities in wake 


1 tj-2 

dynamic pressure, ^pv 


^t 

R 

R<- 


dynamic pressure at tail 


resultant force 


nondimensional radius to origin, 


r*w 


z \2 

+ 


5 


°H 

S V 

u,v,w 


wing area, or equivalent area for nondimensionalizing forces 
and moments 

wing semispan 

sin x H 

sin Xy 

induced velocities along the X-, Y-, and Z-axes, respectively 


u o> v o> w o mean, or momentum theory, values of induced velocities at the 
lifting system 

V forward velocity 


V R 

wh 

X,Y,Z 


x,y,z 




x i’ x j» x kl 

J 


resultant velocity 

reference velocity, vertical induced velocity which lifting system 
would have if it could hover with momentum area Am ’ Wh " ' 

Cartesian coordinates, X-axis positive rearward, Y-axis positive to right 
when viewed from behind, Z-axis positive upward. Unless otherwise 
noted, the origin is centered in the lifting system of the model. 

location of a point with respect to the X-, Y-, and Z-axes, respectively 

location of a point on the wake with respect to the X-, Y-, and Z-axes, 
respectively 

generalized coordinates where i, j, and k may independently take on the 
values x,y,z (Note, for example, that x x = x, x y = y, and x z = z.) 

lateral force 


a 


angle of attack, positive nose upward 


sideslip angle 
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il 


B 

y ratio of wind-tunnel width to wind-tunnel height, — 

6 interference factor (in general terms) 

64 a interference factor for finding the interference velocity in the j -direction 

™ Am 

caused by forces in the i-direction, defined by A^. = 5^ 

®i,j ,1c interference factor for finding the rate of change per semiheight in the 

k-direction of the interference velocity in the j-direction caused by 

f A m 

forces in the i-direction, defined by A ' S i,j,k 

? ratio of wind-tunnel semiheight to height of model lifting system above 

the floor of the wind tunnel 

rj ratio of distance between origin of wake and right-hand (as viewed from behind) 

sidewall of the wind tunnel to the semi width of the wind tunnel, —• 

0 rate of pitch, positive nose upward, rad/sec 


A wing sweep angle, measured positive rearward from lateral axis of 

aircraft, deg 

fjL helicopter tip-speed ratio, ^ C ^ e ed 

v\ generalized mean induced velocities at lifting system, u 0 when i = x, 

v 0 when i = y, w Q when i = z 

p mass density of air or other test medium in the wind tunnel 

a ratio of wingspan to tunnel width, —■ 

13 

$4 a function related to induced velocity in the j-direction caused by forces 

J 

in the i-direction 


. function related to the rate of change in the k-direction of the induced 

velocity in the j-direction caused by forces in the i-direction 



<£= i for a mirror-image wake 

1 jJ 


I 


7 



$i,j,k ^i,j,k f or a mirror-image wake 

<Pi potential of wake as produced by forces in the i-direction 

(Pi j induced velocity in the j-direction as caused by forces in the i-direction 

<?i j rate of change in the k-direction of the induced velocity in the j-direction 

as caused by forces in the i-direction 

Xjj horizontal wake skew angle, angle measured positive rearward from the 

negative Y-axis to the projection of the wake on the X-Y plane 

Xy vertical wake skew angle, angle measured positive rearward from the 

negative Z-axis to the projection of the wake on the X-Y plane 

^ yaw angle, measured positive nose right when viewed from behind 

Prefix: 

A change in value caused by boundary interference 

Subscripts : 

c value corrected for boundary interference 

M value from momentum theory 

i,j, or k indices which may independently take on the value x, y, or z, denoting the 
major axis parallel to which a dimension, force, velocity, or velocity 
gradient is to measured. 


THEORY 

Momentum Considerations 

Certain quantities, upon which the interferences will be found to depend, may be 
found from momentum considerations. The present analysis parallels that of reference 9. 

Consider a force-producing system, acting upon the fluid flowing through an area 
A m and producing lift, induced drag, and Y- force components, as sketched in figure 1. 
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Since force is equal to the time rate of change of momentum 


L = pA m V R (-nw D ) 
Di = pA m V R (-nu 0 ) 
Y = pA m V R (-nv D ) 


where 


V R = ^+Uof+^o 2 + w 0 2 

Dividing equations (2) and (3) by equation (1) yields 

% = ^o 

L w o 

Y = Vo 
L w 0 

Expanding equation (4) and substituting equations (5) and (6) into it yields 



Dividing both sides by -w G 



From equation (1) 


w 0 = 


-L 

npA m V R 


, 2 = -Lw p 
0 npA m V R 




Now define wjj as the induced velocity w Q when V = 0, Dj = 0, and Y = 0; that is, 
w^ is the value of w 0 in a purely hovering condition. From equation (1) 



Divide equation (9) by the square of equation (10) to yield 



npAm vSo) 

L 

npA m 


-w 0 

Vr 


or 



1 



( 10 ) 


(ID 


( 12 ) 


which is best solved for w 0 /wh on digital computing machinery. The desired root is 
generally the smallest positive real root. 

The vertical skew angle (Xy) M > from figure 1, is seen to be 


. . V + u n v 

tan (*V) M = -w Q = ' 


L 



where v/w 0 may be obtained from Wo/wh by means of the identity 

V = V / W h 
w 0 _ w 0 /w h 


The horizontal skew angle 



is similarly seen from figure 1 to be 


(13) 


(14) 


( 15 ) 
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lan Mm 


v + Uo 
_v o 


V_ Uo = w 0 /y u Q \ 
v o~ v o “v 0 \ w o w o/ 


V_ «o x +5i 

w Q w 0 _ wp l_ 

Vq Y 

w o L 


(16) 



(17) 


or substituting equation (13) into equation (16) 


tan Mm 


tan M m 

Y/L 


(18) 


Observe from figure 1 that the wake passes downstream with no deflection whatever 
(as in classical theory) if Xjj and Xy are both 90°. Note that the hovering case, where 
the wake must pass directly downward, is given by Xy = 0° and Xjj = 90°. 


Effective Skew Angles 

Reference 10 observes that the inclinations of the mass flow and of the vorticity in 
the wake differ because of wake rollup. From experimental observations of the wake of a 
rotor (ref. 15), it was concluded that the wake vorticity would be deflected downward by 
approximately one-half the deflection indicated by equations (17) and (18). Under this 
assumption, the effective wake skew angles would be given by 


Xy - 



M),, + 90 ° 

X H 2 


(19) 

( 20 ) 


Equations (19) and (20) suffer from the obvious deficiency that the wake does not 
assume the proper skew angles in hovering. Reference 11 presents an alternate view- 
point which is more aesthetically pleasing in this regard. However, the practical limita- 
tions on very low speed testing (refs. 2, 8, 10, 11, 12, and 14) effectively limit tests to 
wake angles above which there would be any significant differences between equations (19) 
and (20) and the alternate form given in reference 11. The wake angles used in the fol- 
lowing analysis are completely arbitrary and should be interpreted as being the effective 
wake angles, however obtained. 
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Wake in Free Air 


The wake is assumed to be a straight line starting at the model and extending to 
infinity. It is skewed rearward from the vertical axis by an angle Xy and skewed rear- 
ward from the horizontal axis by an angle Xjj (fig- 2). The wake is considered to con- 
sist of a string of point doublets whose axes make some constant arbitrary angles with 
the X-, Y-, and Z-axes. 

Since, in the present linearized theory, the effect of the walls on the wake shape 
is neglected, it is permissible to choose doublets whose axes are parallel to the nega- 
tive X-, Y-, and Z-axes only. Any arbitrary doublet inclination can then be obtained 
as a linear superposition of the above three cases. 

The potential of a single doublet (fig. 2) of the infinite string comprising the wake 
caused by forces in the Xj-direction has its axis directed parallel to the negative 
Xj-axis. Thus, with respect to Cartesian coordinates centered in the point doublet, the 
potential is 


rim . y- 



i 


( 21 ) 


The coordinates of the wake may be expressed as a set of parametric equations 
in l, the length along the wake, as 



1 - cos^Xh cos Xy 


( 22 ) 


-sin x v cos x H 

y = . — - =i 

\/l - cos 2 x h COS 2 Xy 

: -sin x H COS Xy 

z — - J - r- ===== ■ . : l 

/l “ COS 2 X H COS 2 Xy 

More generally, equations (22) to (24) may be rewritten as 


(23) 


(24) 


h = a ± l (25) 

where the ai are defined implicitly by equations (22) to (24). It will be convenient in the 
12 


I 


following derivation to note that ^ ( a i / = * si nce along the wake = x 2 + y 2 + z 2 . 

i 

Now if the coordinate system of the entire wake is chosen to coincide with the origin of 
the wake, the wake potentials may be obtained by an integration over the entire length of 
the wake; that is 



where the summations are considered to be carried out over the three values of i. 

Equation (26) may be integrated immediately by the use of items 162 and 170 of 
reference 16 to yield 



After substituting limits and performing some algebraic simplification, equa- 
tion (27) becomes 



Equation (28) can be written in still more general form as 

dA 

dim 9xi 

*i a ~dTT 

where 

A = lx 2 + y 2 + z 2 ) (l - cos 2 x h cos 2 x v ) 


- x sin x H sin x v + y sin Xy cos X H + z sin x H cos x v 


(29a) 


(29b) 
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When the wake is not deflected sidewards (that is, when Xjj = 90° in fig. 2), equa- 
tion (29) reduces identically to equation (7) of reference 5 for Xj = x, and to equation (1) 
of reference 5 for xi = z. Thus, the theory presented in reference 5 will appear as a sub- 
set of the theory developed herein. 


The induced velocity in the j-direction is then the partial derivative of equation (29) 
with respect to Xj, or 

92 A 

b ( p i 


= n, “ 


dmj 

W 


3XjL 9Xj 


3A 


9A 

9x, 


A2 


(30) 


Similarly, the rate of change in the k-direction of the induced velocity in the 
j-direction is given by the partial derivative of equation (30) with respect to x^, or 


Xi 


dmj 

“d T 


3 2 A 


3x t 8xj 


3A 

0Xi 


3 2 A 


3xr yOxj 3xk J Oxj \dx^ 9x k 

A2 


3A/ 3^ A 


8A / 3 2 A 

Sxjj y3x^ 3xj 


+ 2 


3A 

9X; 


3A \ ( 3A 


3xj / 1 9x, 
A? 


(31) 


Now the values of 


dm^ 

"dT 


are obtained, following reference 5, as 


d m i _ A m 
d l 2ir 1 


(32) 


Further, it will, in general, be convenient to nondimensionalize the field points 
which appear in equations (30) and (31) with respect to h = ^. Substituting equation (32) 
into equations (30) and (31) and performing the indicated nondimensionalization yields 
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where 







( 35 ) 



and 


A 

h 



cos 2 x 


H 


COS 2 Xyl 



) sin x v cosx H 



sin x H cos x v 


(37) 


For convenience, the partial derivatives of 


A 

h 


are given in appendix A. 


It will be desirable to have the appropriate velocities and slopes for a wake whose position (but not strength) is 
a mirror image across the X-Z plane of the original wake. Substitution of 180° - Xjj for Xjj in equations (35) 
►jj to (37) yields immediately 




( 38 ) 


and 



( 39 ) 


where 


F " J j| I ) 2 h) 2 + 1 1) 2 ]^ " cos2 *H cos2 X V ) - (? I) sin X H sin Xy - 1 sin x v cos x H + 1 sin x H cos X y 


( 40 ) 


Those partial derivatives of — which differ from the corresponding partial derivatives of ^ are given in 
appendix B. 



Wake in the Wind Tunnel 


Consider the model located at the center of the coordinate system as shown in fig- 
ure 3 with its wake inclined at angles Xjj and Xy ( a s in fig. 2). 

The wake will intersect the horizontal plane of the floor at a value 


x_ 

H 


z=-h 




Xy 


The wake will intersect the vertical plane of the left-hand wall at 


(41) 


5 ly=-(2B-b) " *“ X H = V8 - 1) % <42) 

Within the physical confines of the test section, the wake will intersect either the 
floor or the wall first. It will be convenient to consider two separate cases: case I, where 
the wake strikes the floor first; and case II, where the wake strikes the wall first. 

Case I: wake strikes floor first .- This case may be distinguished by the fact that 


x 

IT 


z=-h 



[y=- (2B-b) 


(43) 


or 


tan x H > 


tan x v 
?r(2 - v) 


(44) 


In this case the wake descends from the model with skew angles Xjj and Xy until 
it intercepts the floor. At this point it can no longer pass downward and is assumed to 
continue along the floor with the original horizontal skew angle of x H bu t with x v = 90 °- 
Some distance farther rearward, the wake then intercepts the wall at the corner. From 
this point on, the wake passes to infinity in the corner with Xjj = Xy = "°- 

Case II: wake strikes wall first.- This case may be distinguished by the fact that 


x 

H 




z=-h H 


y=-(2B-b) 


(45) 


or 


tan x H < 


tan Xy 
Sr (2 - rj) 


(46) 
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In this case, the wake descends from the model with skew angles x H and Xy until it strikes the wall. At 
this point it can no longer pass sideward and is assumed to continue along the wall with the original vertical skew 
angle of Xy but with Xjj = 90°. Some distance farther rearward, the wake then intercepts the floor at the corner. 
From this point on, the wake passes to infinity in the corner with X„ = X V = 900. 

Interference in Case I: Wake Strikes Floor First 

The wind-tunnel interference in case I is most easily approached by considering two simple systems before con- 
sidering the more involved case of the wind tunnel itself. 

The first simple system is that of a horizontal floor only (fig. 4). This simple system is, of course, that which 
applies to ground effect and will be discussed again in a later section as such. For this case, by superposition from 
the previous results, with the origin centered at the start of the upper wake 



where the following notation has been used 




and where, following in the same line as reference 5, in order to maintain the proper direction of the image strengths 
and image- induced velocities with respect to the coordinate system of the real wake 

p = 1 for i = y, and p = 0 for i = x or i = z 
q = 1 for j = y or j = z, and q = 0 for j = x 


The slopes of the induced velocities for this system may be written as 
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where r = 1 when k = y or k = z, and r = 0 when k = x. 

It is observed that adding slopes by superposition as in equation (48) amounts to adding the tangents of several 
angles to find the tangent of the sum of the angles. Note that this procedure is valid only if <p. . , and all of the indi- 
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vidual terms summed to obtain <p. . . are small. 

Now consider two such pairs of wakes disposed a distance 2B-b to either side of a solid wall as in figure 5. 

For the four wakes, with the origin at the start of the upper right-hand wake, the induced velocity field of the entire 
system is found to be 
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Similarly, the slopes of the induced velocities are given by: 
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The arrangement of images required to insure that there is no flow normal to any of the four walls of a closed 
wind tunnel is shown in figure 6. It will be observed that the repetitive image set which forms the basis for this pat- 
tern is that of the four wakes for which the induced flow field is given by equation (50). This image set is merely 
translated laterally at intervals of 4B = 4y?h and vertically at intervals of 4H = 4£h in order to form the complete 
pattern. 

Therefore, if the interference velocities caused by the presence of the walls are defined as 
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the interference factor 
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The rate of change of the interference velocities follows in the same manner. Defining 
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It will be observed that the central image, which represents the wake in free air, has been omitted in equa- 
tions (52) and (54). This image is removed since it is only the wind-tunnel interferences which are of interest and 
not the total induced fields of model and wind tunnel. 


If Xjj and Xy are both 90°, the wake passes directly rearward without touching the wind-tunnel boundaries. 
Thus, equations (47) to (50), (52), and (54) may be simplified by omitting all terms on the right-hand side of the equa- 
tions except the first, fourth, seventh, and tenth sets of terms. If only x H is 90°, the wake will touch the floor, but 
not the sidewalls of the tunnel. Thus, the third, sixth, ninth, and twelfth sets of terms on the right-hand side of equa- 
tions (49), (50), (52), and (54) may be omitted. If only Xy is 90°, the wake touches one of the walls, but not the 
floor. This set of conditions falls under case n, the equations for which are developed in the following section. 



Interference In Case II: Wake Strikes Wall First 


The wind-tunnel interference in this case is also most easily approached by considering simpler systems first. 
The first such simple system (fig. 7) consists of two wakes disposed a distance 2B - b to either side of a solid wall. 
The induced field of this system with the origin centered at the start of the right-hand wake is 



Similarly, the gradients of the induced velocities are given by 



If the two images are now placed below the preceding two images as in figure 8, the induced velocity field, 
with the origin at the start of the upper right-hand wake, is given by 
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and the gradients of the induced velocities are given by 
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If 6. . is defined as in the preceding section, that is, as 
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It will be observed that the central image describing the field of the wake in free air has been omitted from equa- 
tions (60) and (62). This has been done since it is only the interference due to the walls which is of interest herein. 

If Xjj and Xy are both 90°, the wake passes downstream without touching the wind-tunnel boundaries. Thus, 
equations (55) to (58), (60), and (62) may be simplified by omitting all terms on the right-hand side of the equations 
except the first, fourth, seventh, and tenth sets of terms. (Note that under these conditions in each case, the equa- 
tions are identical in cases I and II.) If only Xy is 90°, the wake will touch the sidewall but not the floor of the wind 
tunnel. Thus, the third, sixth, ninth, and twelfth sets of terms on the right-hand side of equations (53), (58), (60), 
and (62) may be omitted. If only x H is 90°, the wake touches the floor but not the sidewall. This set of conditions 
falls under case I, the equations for which were developed in the preceding section of this paper. 


Interference in Ground Effect 

As mentioned previously, the initial wake image system considered in developing the wall corrections for case I 
corresponds to ground effect. Thus, if the interference in ground effect is expressed as 
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the interference factor can be determined immediately from equation (47) (by taking H = h, so that £ = 1, and y = 1, 
and so that A T = Aq, and by omitting the term corresponding to the real wake in free air) as 
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Note that if Xy is greater than or equal to 90° the wake does not intersect the ground. Under these condi- 
tions, the foregoing equations become simply 
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As in reference 5, it is possible to express the interference factors at the center of 
lift in closed form when dealing with ground effect. Because of the vast multiplicity of 
interference factors involved in the present analysis it is impractical to present complete 
derivations of these factors as was done in the appendixes of reference 5. Instead, only 
the final closed-form expressions are presented in appendix C. The specific closed-form 
expressions for the partial derivatives of ^ and as used in deriving the results of 
appendix C, are given in table I. 


Computer Program 

Program features .- In general, the foregoing equations cannot be evaluated in closed 
form except for a few isolated special cases. Consequently it is necessary to resort to 
high-speed digital computing equipment in order to obtain numerical values. A listing of 
the program used to obtain the values which will be presented in subsequent portions of 
this paper is presented in appendix D. A flow chart for this program is presented in 
appendix E. Other programs were evolved in the course of the study which directly 
obtained values for the average interference, distribution of interference, and interfer- 
ence at the tail for lifting systems which could be represented by arbitrarily swept wings. 
These latter programs are not presented herein since the modifications to the listing of 
appendix D in order to obtain the more involved computer programs will be obvious from 
an examination of references 6 and 7. 

The program of appendix D is very flexible. A series of x H and Xy for which 
values are desired may be inserted as input data. For x H not equal to 90°, an additional 
vertical wake angle is added automatically, this angle being the particular Xy which pro- 
vides the borderline value separating cases I and n. This feature may be eliminated by a 
suitable input character if so desired; it is eliminated automatically if the extra angle is 
within 0.005° of one of the input angles. 

Either ground effect or the closed-tunnel case may be selected by input data. For 
ground- effect calculations, £, rj, and y should be selected as 1.0. 

In addition to the complete solution for all interference factors, two limiting options 
are available. In one, the calculations are limited to only those interference factors which 
pertain to velocities since these nine terms are those of primary interest in simple per- 
formance testing. In the other option, the terms relating to the longitudinal gradients are 
computed in addition to the aforementioned nine terms. This combination provides the 
interference factors of primary interest for tests involving only performance and longitu- 
dinal stability. The selection of either of these two options will reduce the computing 
time substantially. These options are not available for ground-effect calculations since 
the running time is so brief as to be immaterial. 

It is obviously not possible to perform the required summations between infinite 
limits as indicated in the equations of the earlier sections of this paper. After several 
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trial calculations, it was found that summing over m and n between the limits of ±3 
appeared to achieve a reasonable compromise between speed and accuracy. Since the 
actual test section occupies one corner of a repetitive group of four images, this choice 
results in six columns of images to the right and seven columns of images to the left (as 
viewed from behind) of the real test section. Similarly, there will be six rows of images 
above and seven rows of images below the real test section. This image pattern differs 
from that of references 4 to 8 wherein there were only three image sections on each side 
of the test section, but the vertical arrangement therein is identical to that of the present 
paper. 

Numerical checks .- Within the limits imposed by the different image systems, the 
interference factors obtained herein (when x^ = 90°) are identical to those of references 4 
to 8 which, in turn, are completely compatible with the values obtained (xy = Xjj = 90°^ by 
more conventional wall-interference calculations. Small residual values are obtained 
at Xjj = Xy = 90° for certain terms for which symmetry requires that the values be zero. 
This result occurs because of the lack of complete symmetry in the image systems of the 
present paper and also in the image systems of references 4 to 8. These residual values 
could be reduced by increasing the limits of the summations; however, the increased com- 
puting, time should not be worth the minor increase in numerical accuracy. 

Several of the gradients can also be compared with references 5 to 7 by comparing 
the gradient computed herein with the plotted results for a series of points computed by 
the procedures of the earlier papers. Numerous such comparisons have been made, and 
in all cases the slopes agree. 

The foregoing checks on numerical accuracy do not encompass cases in which the 
wake may be deflected to the side as well as downward. There are, however, numerous 
symmetries and equivalences which must be met in the results. These features will be 
discussed in the following several sections of the paper. All of these additional tests are 
met by the present theory as implemented by the program given in appendix D. 

In ground effect it is possible to obtain a numerical check of the program accuracy 
for values of Xjj other than 90° since closed-form solutions for the interference factors 
at the center of lift have already been presented in appendix C. A comparison of the com- 
puted values obtained in both manners indicated complete agreement. 

An additional numerical test was performed by calculating the interferences for a 
small model in the center of a square tunnel and choosing selected equal values of Xjj 
and Xy These choices duplicate the wake and tunnel configuration of the "diamond" test 
section of reference 17, for which results obtained by a different theoretical treatment are 
available. After the resolution of interference vectors required by the differing coordinate 
systems, the present results are identical to those of reference 17. 
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Skew angles in the second quadrant .- It has already been noted that the image sys- 
tems in ground effect differ according to whether or not the vertical wake skew angle is 
greater than 90°. Thus, in ground effect, the program automatically chooses between 
equations (65) and (67) or equations (68) and (69) according to the value of Xy 

Although it may not be as evident, the equations developed for the closed tunnel will 
also lead to an incorrect wake pattern if either x H or Xy is greater than 90°. As will 
be discussed subsequently, appropriate values, for skew angles exceeding 90°, can be 
obtained from certain symmetry relations. This procedure is embodied in the computer 
program, which automatically sets up an equivalent problem in the first quadrant and then 
converts the results to correspond with the required angles in the second quadrant. 

The superposition techniques used herein to obtain the field of wake segments of 
finite lengths can lead to numerical difficulties (of the nature of overflow and underflow) 
if the point of interest in the tunnel lies on or near an extension of any of the image wakes. 
Particular difficulty will be experienced for wake angles near zero. The program as 
presently constituted excludes values of the skew angles which are less than or equal to 
zero. Minor modifications to the program at lines (D 108) and (D 111) will allow the use 
of angles less than but not equal to zero; however, if zero is given as an input value, the 
execution of the program will terminate at that angle. Since the results of references 8, 

10, 13, and 14 indicate that reasonable testing conditions cannot be obtained at very small 
skew angles, the practical effect of this restriction should be minimal. 

Program storage requirements and running time .- The program of appendix D as 
implemented in the CDC 6600 computers of the Langley Research Center requires 42 000g 
(approximately 17 500 jq) spaces in memory in order to compile. Execution of the program 
requires 24 000g (approximately 10 300 iq) spaces. 

The time required for the calculations varies according to the particular grouping of 
coefficients desired, the boundary conditions (ground effect or closed tunnel), and, since it 
is permissible to omit certain terms in those cases, whether Xjj> Xy> or both are 90°. 

The following table gives the approximate central processor time (in seconds) required for 
each combination of one x H and one Xy : 


Approximate central processing time, in seconds, in the Langley computer complex 


Boundary 

conditions 

Coefficients 

General 

X R = 90° or Xy = 90° 

X H - Xy = 90° 


All 

9.3 

5.7 

2.1 

Wind tunnel 

6 i,j and 5 i,j,x 

5.7 

3.5 

1.3 


5y only 

1.8 

1.1 

0.5 

Ground effect 

All 

0.8 

0.4 
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In addition, when the computer task includes compilation from the FORTRAN listing, 
approximately 3.7 seconds will be required for the compilation. 

Program efficiency .- The present computer program is not particularly efficient, 
and the computer time required when the program is expanded (by the techniques of ref. 6) 
to finite-span configurations is clearly excessive for routine application. In a sense, this 
inefficiency was deliberate in that it was desired to compute all terms in order to obtain a 
check on the symmetries rather than to use the symmetries to reduce the required running 
time. Substantial reduction in computer time could be obtained by reprograming the calcu- 
lations to take advantage of the symme tries. 


Symmetry of Interference Factors 

Wake in free air .- Because A and A' are both jointly continuous through at least 
the third order, as can be seen from equations (37) and (40), transposing the i and j in 
<3>i,j (eq. (35)) or 4>j j (eq. (38)), or transposing i, j, and k in any possible permutation 
in 4>ij t (eq. (36)) or ( e< h (39)) leaves these terms unaltered. Thus, there is a 

possibility that certain of the interference factors will also be unaltered by such opera*- 
tions. This portion of this report considers these symmetries separately for the cases of 
the wind tunnel first and then for the case of ground effect. 

Wind-tunnel interference .- Examination of equations (52) and (60) indicates that an 
additional requirement for equality of the correction factors for the interference velocities 
is that (-1)^, (-1)9, and (-lJP+Q must also be unaltered if the correction factors are to 
be unaltered by transposing i and j. No combinations of i and j exist which satisfy 
these requirements; therefore, in general, all of the correction factors for interference 
velocities differ. On the other hand, examination of equations (54) and (62) indicates that 
the additional requirement for equality of the factors for the slopes of the interference 
velocities is that (-l)P, (-l)9 +r , and (-l)P + 9 +r be unaltered. Because of the way in 
which p, q, and r are related to the i-, j-, and k-directions, this requirement is far less 
restrictive than the corresponding requirement for the velocities, and immediately yields 


6 


i,j,k - 6 i,k,j 


(70); 


Ground effect .- In ground effect, conditions for the equality are somewhat less 
restrictive. For example, for the interference factors describing velocities, it is only 
necessary that (-ljP 4 ^ be unaltered if the factor is to be unaltered by interchanging i 
and j. Thus, 


6 x,y 6 y,x 


(71) 


Similarly for the interference factors describing slopes, it is only necessary that 
(-l)P+q+ r be unaltered by a permutation in i, j, and k. Thus, in addition to equation (70) 
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6 =6 =6 
x,y,y y,y,x y,x,y 

(72a) 

6 =6 =6 
x,x,y y,x,x x,y,x 

(72b) 

6 =6 =6 =5 

x,y,z y,x,z x,z,y y,z,x 

(72c) 


The foregoing symmetries are displayed in table II, which reproduces the computer 
output for one set of conditions. 


Interchange Equivalences 

Between xh an d XV- " Ah interesting and sometimes useful set of equivalences may 
be found by considering the manipulations shown schematically in figure 9. Consider an 
arbitrary initial wind tunnel as in figure 9(a) with the wake deflected both laterally and ver- 
tically, and at arbitrary point (x,y,z) at which the interferences are known. Rotation of the 
entire picture 90° counterclockwise yields the configuration shown in figure 9(b). Then 
reverse the configuration to a mirror image as in figure 9(c), and, finally, relabel the con- 
figuration to be in accord with the standard symbol nomenclature of this paper. The fol- 
lowing table lists the pertinent quantities relating to figure 9(d) in terms of the initial val- 
ues of figure 9(a). 


Initial tunnel 

Derived tunnel 

X H 

Xy 

Xy 

X H 

x/H 

(l/y)(x/H) 

y/H 

(l/y)(z/H) 

z/H 

(i/y)(y/H) 

y 

l/y 

S 

1/(2 - V) 

V 

2 - I/? 


Since it is obvious that there is no essential difference in the physical tunnels if the 
above conditions are met, the interference factors will be identical in the two cases pro- 
vided only that in the subscripts of , z and y are interchanged for each other wher- 
ever they occur. For example, in the final tunnel (fig. 9(d)), 6 X>Z , 6y,y, 311(1 6 z,y are 
identical to 6 X> y , 6 z ,z> and 5y jZ , respectively, in the initial tunnel. 

The gradients are only slightly more involved. In terms of constant unit lengths 
these quantities would be related in precisely the same manner as the velocities; however, 
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the gradients as obtained herein are with reference to the nondimensional quantities Xj/H. 
Thus, regardless of the direction of the gradient, its dimensions are those of a velocity per 
unit semiheight. Since the semiheight of the original tunnel is interchanged with the semi- 
width of the final tunnel, these factors may be corrected to the new semiheight by multi- 
plying the interchanged Si^k by y (that is, for example, 6y,x,z in the final tunnel is 
equal to (y) • (6 z ,x,y) * n the initial tunnel). 

These relationships are demonstrated in table in by direct reproduction of the com- 
puter output for cases meeting the foregoing criteria. These relationships are particu- 
larly obvious for a centrally located model in a square test section, and still more obvious 
when Xjj and Xy are identical. 

The foregoing interchange relationships were very useful in checking the program of 
appendix D for cases in which the wake is deflected laterally. Not only do they provide a 
check on consistency of the results, but for Xy = 90° the computed interferences may be 
compared directly with the computed values of reference 5 for the equivalent tunnel. 

It is observed that for wind-tunnel stability tests of conventional aircraft configura- 
tions, where the wake may generally be considered to be undeflected in either direction, 
these interchange equivalences may allow one to obtain the correct interference factors 
for side forces directly from available information (ref. 1) for 'lift interference.” 

Strictly speaking, the interchange equivalences apply only to vanishing small 
models. If the model had been considered as a finite-span wing lying on the Y-axis of the 
initial tunnel, figure 9 indicates that the span would be lying on the Z-axis of the final tun- 
nel. The net result would be an entirely different spatial distribution of vorticity in the 
tunnel. Under such circumstances, there is no reason to expect that the interference fac- 
tors would follow the rules given here. Several sample calculations indicate that the theo- 
rem is not greatly violated for relatively short spans (say ct = 0.25); however, for large 
spans, significant differences arise (table IV). 

Between values of xh in th e first and second quadrants.- Now consider figures 9(e) 
and 9(f). Except for the positive direction of the Y-axis in figure 9(f), these two figures 
are equivalent provided that the conditions given in the following table are met: 


Figure 9(e) 

Figure 9(f) 

X H 

180« _ XH 

7 

2-7? 

y 

-y 


The effect of changing the positive sense of the Y-axis with respect to the wake 
depends solely on the number of differentiations with respect to y/H that are implied by 
the subscripts of ^i ? j, $i,j,k> 6 i,j> and Si j k since, for complete equality of the results 
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between the two figures, it would be necessary to differentiate with respect to -y/H for 
figure 9(e). Thus, the interference factors in which the subscript y appears an even 
number of times will be identical in the two cases, and those in which y appears an odd 
number of times will be identical except for a factor of -1. As discussed earlier, this 
interchange is utilized in the FORTRAN program of appendix D to obtain interference fac- 
tors when x H lies in the second quadrant. 

If the model is in simple ground effect, or if the model is centered laterally in the 
closed tunnel, this interchange requires that the interference factors in which the sub- 
script y appears an even number of times be symmetrical with respect to Xjj = 90°. 
Thus, this group of interference factors will be referred to as the symmetric factors 
herein. Similarly, if y appears an odd number of times, the interference factors will 
be antisymmetric about Xjj = 90°; consequently, this group will be referred to as the 
antisymmetric factors. 

Observe that for laterally centered models a symmetric factor may have a substan- 
tial value at Xjj = 90°, but that it will vary relatively slowly for small departures from 
90° because the rate of change of the factor must be zero at that point in order to allow 
symmetry. Since most lateral-directional stability testing involves only small lateral 
deflections of the wake, it may be perfectly acceptable to ignore the effect of Xu on 
these factors, using only the values obtained at Xjj = 90°. 

Further, in order that antisymmetry be maintained, all of the antisymmetric group 
of factors must be zero at Xjj = 90° if the model is laterally centered. On the other 
hand, the rate of change of any of the factors with respect to Xjj could assume a large 
value. Thus some caution must be used in applying the present results if it is arbitrarily 
decided that the lateral deflection of the wake will be ignored. 

In the special case of a laterally centered model with Xjj = 90°, the conditions spec- 
ified in the foregoing table require that the symmetric factors have a symmetric distribu- 
tion across the Y-axis. Under the same conditions, the antisymmetric factors will have 
an antisymmetric distribution across the Y-axis. 

Between values of Xy * n th e first and second quadrants.- Finally, consider fig- 
ures 9(g) and 9(h). Except for the positive direction of the Z-axis in figure 9(h), these 
two figures are equivalent provided that the following conditions are met: 


Figure 9(g) 

Figure 9(h) 

Xy 

180° - XV 

s 

?/(2? - 1) 

z 

-z 
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The effect of changing the positive sense of the Z-axis is similar to the effects dis- 
cussed in the preceding section. The interference factors fall into symmetric and anti- 
symmetric groups depending upon whether z is repeated an even or an odd number of 
times in the subscripts. Symmetric interference factors will be unaltered by the inter- 
change, and antisymmetric factors will be altered in sign only. 

If the model is centered vertically in the tunnel, the symmetric factors will be 
symmetric about Xy = 90° and the antisymmetric factors will be antisymmetric about 
Xy = 90°. Again the symmetric factors may have significant values at Xy = 90°, but 
these factors will vary only slowly for small deviations from 90°. The antisymmetric fac- 
tors must be zero when Xy = 90°, however, there is a possibility that they may be signifi- 
cantly altered for comparatively small vertical deflections of the wake. Furthermore, 
when Xy = 90°, the symmetric factors will have symmetric distributions over the Z-axis 
and the antisymmetric factors will have antisymmetric distributions over the Z-axis. 

It is important to observe that this interchange is completely voided in ground effect. 
The reason, of course, is that there is no restraint on the flow above the model. The only 
restraint is the ground below the model. 

APPLICATION OF RESULTS 

Wall-interference theory, in general, is limited to the calculations of the interfer- 
ence velocities and gradients which the walls contribute to the overall flow in the wind 
tunnel. Ideally, the process of correcting wind-tunnel data involves calculating the effect 
of these velocities and gradients upon the model and then subtracting these effects to 
obtain interference-free data. This ideal process is seldom attempted in practice. 

Indeed, for some of the more exotic V/STOL vehicles, there is no background theory 
available with which to make such calculations. In other cases, such as wings (for which 
modern vortex- lattice techniques could be used for relatively accurate calculations), the 
required computer time for elaborate correction techniques is excessive, and corrections 
to data are generally made by far more approximate techniques (ref. 1). 

An alternative to the direct calculations and removal of wall-interference effects is 
to consider data correction as a problem in similitude. This is the most generally applied 
technique. In this manner, the performance as measured in the tunnel becomes simply 
the correct performance for an altered flight condition, or even for a model slightly 
altered from that actually tested in the wmd tunnel. At times, such similitudes actually 
lead to alternative means of approximate direct calculations. This portion of the present 
paper will discuss a number of such techniques which, hopefully, will be of some help in 
applying the results of the foregoing analysis. 
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Interference Velocities 


The preceding portions of this paper have obtained the interference velocities along 
each of the three axes of the tunnel as caused by forces along the three axes. The com- 
plete interference velocities are therefore given by 



(73a) 

(73b) 

(73c) 


Altered flight velocities .- The interference velocities given by equations (73) are 
superimposed on the main wind-tunnel velocity. Conceptually, the simplest means of 
accounting for these interference velocities is to assume that the measured data corre- 
spond to the performance in free air of a model having a forward speed of V + Au and a 
sink rate of Aw while it is translating to the left at a velocity of Av. When the main tun- 
nel velocity is small, this technique may be reasonably satisfactory since it corresponds 
roughly to the conditions occuring during a landing approach in a cross wind. At speeds 
corresponding to hovering (or near hovering) conditions it may even be the only reason- 
able approach; however, in general, in order to obtain meaningful equivalent flight condi- 
tions, it will be convenient to express at least one of the interference velocities as an 
angular change at the model. 

Altered yaw angle .- Consider first an equivalent flight condition in which the longitu- 
dinal and vertical interferences are still considered as mere alterations in linear speeds, 
but in which the lateral interference velocity is to be expressed as an alteration of the 
effective yaw angle of the model. From figure 10, it may be seen that the effective for- 
ward speed becomes 

V c = j/(V + Au) 2 + (Av) 2 (74) 

and the corrected dynamic pressure is 

q c =|p|(V + Au) 2 +(Av) 2 ] (75) 

All of the coordinates in the tunnel-based axis system must now be altered to corre- 
spond with the new axis system by means of a rotation through an angle tan - 1 |av/ (V + Au)J 
about the Z-axis. In matrix form 
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x c 


V + Au 

n 


X 


|/(V + Au) 2 + (Av) 2 

/(V + Au) 2 + (Av) 2 



y c 


-Av 

V + Au q 


y 


fcv + Au) 2 + (Av) 2 

/(V + Au) 2 + (Av)^ 


z c 


0 

0 1 


z 


( 76 ) 


The corrected yaw angle is obtained by transforming a unit vector along the forward 
axis of the model into the new coordinate system, and then observing that 


tan^ c = 

c x c 


(77) 


Substitution of equation (76) into equation (77) yields 


tan » p c = 


Av , , 

V + Au + tan ^ 


Ay 


V + Au 


tan rf/ 


or 


where 


- V + A& 


A^ = tan"! LAv 

\V + Au 


(78) 

(79a) 

(79b) 


It is easily verified that a is unaltered by the rotation since the rotation is about the 
Z-axis. 

In terms of the correction angle A^, the matrix transformation given by equa- 
tion (76) becomes 


x c 


yc 

= 

Zc 



cos A 4/ 
-sin A\j/ 
0 


sin A ip 0 
cos A ip 0 
0 1 


(80) 


In the new coordinate system the force coefficients are obtained by subjecting the 
resultant-force vector to the same transformation as the simple coordinates (eq. (80)), 
and then normalizing the force components with respect to the corrected dynamic pressure 
given by equation (76) to obtain 
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C D C = ( C D cos A ^ + Cy sin 

C y c = (-C D sin Axjy + Cy cos A ^)j(^ 


where 


CL c ’ Wit 


T=( 1 + f ) 2 + ff ) 2 


In a similar manner, the components of the corrected moment are obtained as 


C l c = ( c l cos A ^ + c m sin A*)/f 


Cm c = (-Cl sin A^ + C m cos A ^)/ bf 


(82c) 


Altered yaw angle and angle of attack .- An alternate viewpoint is to consider that the 
effective axis of the tunnel airstream is altered by the wall interference and that an equiv- 
alent flight condition may be obtained by merely resolving all of the data about the cor- 
rected stream axis. This is, by far, the most usual manner of dealing with wall effects. 
Referring to figure 11, it is obvious that the effective forward velocity along the new axis 
system will be 


V„ = i /(V + Au) 2 + (Av) 2 + (Aw) 2 


or, in terms of dynamic pressure 


1 +^ + 


All of the coordinates in the old tunnel-based coordinate system must now be altered 
to correspond with the new effective stream axes. The transformation may be 
accomplished by rotating the axis system about the Y-axis (fig. 11) through an angle 

tan -1 [Aw/(V + Au)] and then about the new vertical axis through an angle 

tan -1 A yj fy + Au) 2 + (Aw) 2 J. Substantial complications are contributed by the 
fact that only one of the foregoing axis transformations may be taken about the 
original wind-tunnel axis system. 
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In matrix form 



or, after multiplication of the transformation matrices, 



£» 

00 


( 86 ) 



Expansion of equation (86) into its scalar components yields 


(V + Au) X + Avy + Awz 


/(V + Au)^ + (Av)^ + ( Aw)^ 

(87) 

-Av(V + Au)X + [(V + Au) 2 + Aw^jy - (Av Aw)z 

(88) 

/(V + Au) 2 + (Av) 2 + (Aw) 2 /(V + Au) 2 + (Aw) 2 


-(Aw) X + (V + Au)z 

(89) 

+ Au)^ + (Aw)^ 



The corrected angle of attack a c and yaw angle i^ c (=-/3 c ) are obtained by trans- 
forming a unit vector (-cos a cos ip)i + (cos a sin I//) j + (sin a?)k along the for- 
ward longitudinal model axis into the corrected axis system, and then noting that for 
this unit vector 


tan o>c = 



z 


c 



2 


and 

tan ^ 

Performing the indicated operations on equation (90) yields 


(90) 


(91) 


tan a c 


Aw cos a cos if; + (V + Au) sin a 



(92) 


which may be simplified by using the relationships of the right triangle to obtain 


sin oi c = 


_ Aw cos a cos + (V + Au) sin a 
y/(V + Au)^ + (Aw)^ 


When p = 0, equation (92) may be reduced to 


(93) 
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(94) 


tan a + 


Aw 


tan a Q = 


1 - 


Aw 

V + Au 


— + - Au = tan [a + tan -1 


tan a 


Aw 


V + Au 


Therefore, when ^ = 0 


a c = a + Aa 


where 


A a = tan' 


■1 Aw 


(95a) 


(95b) 


V + Au 

Similarly, performing the aforementioned operations on equation (91) yields 

Av(V + Au) cos oi cos ip + [(V + Au) 2 + Aw^cos a sin - Av Aw sin a 


tan = 


|/(V + Au) 2 + (Aw) 2 [(V + Au) cos a cos - Av cos a sin - Aw sin a] 

(96) 


For the special case where a - Aw = 0, equation (96) reduces to 


ta n * = X3M + ^ t = tan L + tan-1 

' V+Au 


(97) 


or 


xfy c = xjy + A\fy 


(98a) 


where 


A\p = tan 1 — — — (98b) 

V + Au 

Next, the force coefficients must be resolved about the corrected axis system. 

This may be accomplished by transforming the resultant-force coefficient 
Cr = i Cr + j Cy + kCR in a manner identical to that of equation (86) and non- 
dimensionalizing the result with respect to the corrected dynamic pressure 

q c = ^p|(V + Au ) 2 + (Av ) 2 + (Aw) 2 J; thus, 
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(99) 


Ct> - — 


V 2 j(V + Au)C d + Av C y + Aw Cl] 
[v + Au) 2 + (Av) 2 + (Aw) 2 ] 


3/2 


V 2 (-Av(V -t- Au) Cd + [(V + Au) 2 + (Aw)^]cy *- Av Aw Cl) 

Cy = - — - - 

C |v + Au) 2 + (Av) 2 + (Aw) 2) 3/2 /(V + Au) 2 + (Aw) 2 


( 100 ) 


C L c = 


V 2 EawCd + (V + Au)cJ 


(V + Au) 2 + (Av) 2 + (Aw) 2 ] /(V + Au) 2 + (Aw) 2 


( 101 ) 


Observe that the lift-drag ratio is also altered; that is, divide equation (101) by equa- 
tion (99) to obtain 


±\ - 

-Aw Cd + (V + Au) Cl 

/(V + Au) 2 + (Av) 2 + (Aw) 2 

l 

(V + Au) Cd + Av Cy + Aw Cl 

J (V + Au) 2 + (Aw) 2 


(102) 


In addition, the moments must also be resolved about the corrected axis system. 
Following the same procedure as for the forces, rotate the resultant moment vector into 
the new coordinate system and nondimensionalize with respect to q c to obtain 


_ V 2 [(V + Au) Ci + Av C m + Aw C n ] 
[(V + Au) 2 + (Av) 2 + (Aw) 22 ^ 2 


V 2 { 

- Av(V + Au) Ci + [(V + Au) 2 + (Aw)^] 

C m - Av Aw C n j 


[or + Au) 2 + (Av) 2 + (Aw)^j 

/ /(V + Au) 2 + (Aw)^ 


(103) 


(104) 


C n c " 


V 2 [ Aw Ci + (V + Au) C n ] 


J]v + Au) 2 + (Av) 2 + (Aw) 2 ] /(V + Au) 2 + (Aw) 2 


(105) 


Scalar qualities .- Scalar quantities such as power are unaltered by the foregoing axis 
transformations; however, if such quantities are formed into nondimensional terms by 
dividing by V or q, then the nondimensional parameter must be corrected to correspond 
to the corrected V or q. Thus, for example, if 
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(106a) 


_ Momentum 
^ qS 


then 



cyv 2 

[(V + Au ) 2 + (Av)^ + (Aw)^j 


(106b) 


Particular caution is required in correcting certain other parameters; for example, the 
tip-speed ratio used in helicopter tests 


_ V cos a 
^ Tip speed 


(107) 


requires corrections to both velocity (eq. (83)) and angle of attack (eq. (92) or (93)). 


Interference Gradients 

General comments. - If the interference velocities were completely uniform over the 
model, there would be little difficulty in correcting wind-tunnel data; indeed, the foregoing 
several sections of this paper would suffice for an almost complete treatment. If the 
model was extraordinarily small, the deviation of the interferences from the average val- 
ues would be negligible, and thus, in combination with the small ratio of model to tunnel 
sizes, it would be permissible to neglect the effect of the interference velocity differences 
over the model. Unfortunately, in the usual wind-tunnel tests, the model has significant 
dimensions compared to the wind-tunnel test section, and significant variations in wall- 
induced interference velocities occur over the extent of the model. Thus, it is necessary 
to account for the effect of these differences in interference velocity over the model. 

The actual techniques by which corrections for interference nonuniformity could be 
calculated depend entirely upon the available theoretical aerodynamic treatments in non- 
uniform flow for the particular model under test. For a conventional airplane configura- 
tion, the powerful numerical techniques of vortex- lattice theory, at least in principle, 
could be used to estimate the effects of interference nonuniformity and to develop correc- 
tions to the measured data. Unfortunately, this approach might possibly involve such 
excessive amounts of computer time that it would be uneconomic for routine use, and even 
then would not account for the possibility that the nonuniform wall-induced interference 
might alter flow separations occurring as the model conditions approach either complete 
or local stall. For less conventional models, this technique is not generally available, 
and for many of the profusion of V/STOL aircraft types, there is no truly adequate theory 
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at all. Thus, in most cases, it will be necessary to use correction treatments which are 
only crude approximations of the actual effect of flow nonuniformity. 

Even if the most rigorous type of theoretical treatment were available for correcting 
data, a certain amount of caution would be required. If the size of the corrections becomes 
comparable to the magnitude of the measured data, it is obvious that the corrected data are 
as much a product of theory as of the wind-tunnel measurements. At this point it becomes 
questionable as to whether or not the wind-tunnel tests have contributed significantly to the 
determination of the model characteristics. 

The surest procedure which avoids the problems inherent in large corrections for 
flow nonuniformity is to avoid testing models which are large relative to the test-section 
dimensions. Many rules of thumb have been offered as to allowable model sizes; however, 
none of these rules is failure free. This trend is enforced by recent studies which indi- 
cate limits to theoretical calculations because of wake deformation (ref. 18) and the need 
for a reduction in allowable correction size when testing swept wings (ref. 19). When the 
primary concern is nonuniformity, a reduction in model size is doubly helpful: first, the 
reduction in area ratio reduces the overall corrections immediately; second, the model 
does not extend as far in any direction and thus, for a given rate of change of interference 
in that direction, will experience a lesser nonuniformity. 

If the model size is reasonably small, many of the interference gradients are rela- 
tively uniform across the model. Thus, some simple concepts can be used to develop com- 
paratively simple first-order corrections. Thus, it is profitable to discuss a few of these 
concepts in the following pages. For the sake of convenience in the discussion, the model 
axes will be considered to be coincident with the principal axes of the tunnel. Similar con- 
siderations would apply for any other model orientation, although considerable resolution 
of the interference vectors might be required to obtain an equivalent frame of reference. 

The direct gradients .- Consider first the rates of change along each axis of the inter- 
ference velocities directed along those axes; that is, 


9(Au) _ y A 

3 (x/H) ” L ^’ x ’ x 

i 

3(Av) _Y a 
3(y/H) 2 A<p i,y,y 
i 


(108a) 


(108b) 



i 


(108c) 


A positive value of any of these three derivatives indicates that the interference 
velocity increases with positive distance along the axis. Compared to a uniform stream 
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in free flight the equivalent free-air model would have a velocity field which would expand 
along these axes. While such a concept seems difficult, it is precisely equivalent in the 
X-direction to the q^/q term applied to tail forces in most stability analyses. For con- 
ventional models, it would be expected that the largest effect of these derivatives would be 
an alteration in the effective dynamic pressure at the tail. 

The gradients 9( Aw)/ 9(x /H) and 9(Au)/9(z/H) .~ Consider next the gradients 



(109a) 


9 !> u ) = y A( p . 

9(z/H) L 
i 


(109b) 


The positive sense of these gradients is indicated in figure 12(a). Along the X-axis 
the distortions in the vertical components of the interference velocities are the same as if 
the model in free air were rotating about the Y-axis in a nose-up direction. On the other 
hand, the horizontal velocities along the Z-axis are those that a model in free air would 
experience if it were rotating about the Y-axis in a nose- down direction. Furthermore, in 
consequence of the symmetries of the interference factors = 6^^, the two equiva- 

lent rotational rates must be identical except for their opposing directions. 

If the model was spherically symmetrical, and if there was no forward velocity, such 
a distortion of the flow field might have little or no effect on the observed performance. 
Such conditions are seldom met, however, and the effects encountered in practice may be 
approximated by considering separately the different portions of the model. If the model 
was a conventional aircraft and lay in the X-Y plane, the effect would be primarily that of 
an effective pitch rate at constant angle of attack; that is (from ref. 2), 


§ - L A(Aw) 
H 9( X /H) 


( 110 ) 


The concept of an equivalent rotational rate may be helpful in developing corrections 
since a large number of studies exist in aerodynamic stability theory for the estimation of 
the effects of such rates of rotation. (See, for example, the survey presented in ref. 20.) 
An alternate viewpoint might be to try to account for the combined effects of an altered 
wing camber, tail height, and tail incidence (ref. 2) without any rotation rate. 

If the model had a significant vertical extent (for example, a biplane of large gap), 
it might well be necessary to consider separately the effects of the tail system rotating 
at 6 and a biplane cell rotating at - 6. Further, if the model had a high T-tail, it might 
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be necessary to consider the additional alteration in the effective longitudinal velocity at 
the tail which would result from the vertical displacement above the plane of the lifting 
system. 

For some simple systems, a search of older studies may even produce means of 
calculating the effects of these gradients directly. For example, reference 21 examines 
directly the effects of a linear gradient of induced velocity normal to the plane of a cen- 
trally hinged rotor and finds that the effect is largely that of an altered lateral flapping. 
Observe that if such a rotor is tested at extreme negative angle of attack the direction of 
the effects will change as the rotor plane approaches the vertical. It should also be 
observed that if the rotor hinges are eliminated and the blades are very rigid, that the 
effect of this gradient changes from a relatively innocuous lateral tilt to a very powerful 
pitching moment (ref. 22). 

It is important to observe that entirely different corrections are appropriate to dif- 
ferent configurations. For some models many of the interference velocities and gradients 
calculated herein will have no significant effect whatever and thus can be ignored. This 
does not necessarily mean that they can always be ignored safely since comparatively sim- 
ple changes in configuration may affect the relative magnitudes of the different effects. 

The gradients 9(Au)/9(y/H) and 9(Ay)/9(x/H) .~ The gradients 





(111a) 


and 


8 (Ay) _ y 
•8(x/H) L 


<p- 


i,y,x 


(111b) 


are shown schematically in their positive senses in figure 12(b). Again the two effective 
rotations are equal but in opposite directions because of symmetry. 

If the model was a conventional aircraft, the effects would be much the same as if 
the wings were experiencing a constant rate of yaw while the fuselage and tail would be 
experiencing the same rate of yaw in the opposite direction. The overall effect on the 
yawing moment corrections would depend upon the balance of these two effects. 

If the model is a helicopter rotor, the effective rotational speed is increased when a 
blade lies along one axis and decreased when it lies along the other. There is no change 
in the average effective rotational speed of the rotor even though a small ripple is imposed 
on this rotational speed. 
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The gradients 9(Ay)/9(z/H) and 9 (Aw) /9(y/H) .~ These gradients, shown schemat- 
ically in figure 12(c) are given by 


9(Av) V 

9^7h) ~Z n >y> z 

i 



(112a) 


(112b) 


Again because of the symmetry of the interferences the two gradients are equal but 
opposite in direction. Thus, an axisymmetric model, such as a propeller, located in this 
plane would have no effective change in rotational rate. On the other hand, a wing lying 
along the Y-axis would be subjected to a wall-induced flow which would be the equivalent 
of a constant rate of roll in free air. 

One aspect of dealing with distortions on the basis of a simple gradient should be 
noted carefully. Observe that for the wing on the Y-axis, only the antisymmetric part of 
the distortion can be described by a single gradient at the center of lift. Symmetrical dis- 
tortions are entirely lost. Thus, if the interference increases or decreases symmetrically 
outward from the center of lift (as is usually the case with at least part of the interference 
in the lateral direction) this information will be lost. Such symmetric vertical interfer- 
ence along the Y-axis is equivalent to an effective alteration (with respect to free air) of 
the wash- in or wash-out of a wing (ref. 2). Similarly, a lateral gradient of 9(Aw)/9(y/H) 
would be equivalent to an effective alteration of the spanwise distribution of airfoil- section 
camber. Such wall-induced distortions can have significant effects on several observed 
quantities (such as the stall angle) and if the wing is swept can have significant effects on 
the longitudinal pitching moment as well (ref. 6). In correlation studies in which tests of 
the same model in different tunnels are compared (such as ref. 23) the pitching moments 
often correlate less well than the lifts and drags. The reason for the poorer correlation 
of moment can often be found in the failure to account for the effect of wall interference on 
the spanwise load distribution of a swept wing or a failure to account for the differences in 
wall interference at the wing and the tail (as in ref. 24). 

If the model has any significant size with respect to the tunnel dimensions, it will be 
necessary to account for this size both by superposition of the present results to obtain 
interference factors corresponding to the actual finite-size configuration and by using the 
actual distributions of interference over the model in order to develop proper corrections 
from the interference velocity calculations. Some aspects of this problem will be exam- 
ined in a later portion of this paper. 
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NUMERICAL RESULTS 


Since the interference factors are expressed herein in an open form, it is difficult to 
determine the effect of changes in the input variables by examination of the equations them- 
selves. Consequently, numerical results obtained from the previously described computer 
program will be presented to illustrate certain effects on the interference factors. 

The only wind-tunnel configuration for which the present numerical results were 
obtained has a width- height ratio y of 1.5. This width-height ratio approximates that of 
many operational wind tunnels. The interchange symmetries discussed earlier allow an 
immediate conversion of these values to correspond to a tunnel having a width-height ratio 
of 2/3. The model is centered in the tunnel at all times. 

Interference at Model 

Effect of Xy" Figure 13 compares the interference factors in the wind tunnel and 
in ground effect for a wake which is laterally undeflected (Xjj = 90°). As predicted by 
the prior consideration of the effect of interchanging horizontal skew angles between the 
first and second quadrants, all of the interference factors for which the subscript y is 
repeated an odd number of times (the antisymmetric group with respect to y) are zero 
throughout the entire range of vertical skew angles. The remaining group (or the sym- 
metric group with respect to y) displays a pronounced dependence upon the vertical skew 
angle with large values being obtained at low vertical skew angles. 

The trends of the symmetric group of factors are similar in both the wind tunnel and 
in ground effect. Note that the factors for ground effect are defined with respect to an area 
(Aq = 4h 3 ) which is only two-thirds the reference area = 4BH = 4yH 3 ) with respect to 
which the factors for the wind tunnel are defined. For small vertical skew angles, this 
factor is the largest part of the difference between the corresponding factors in the tunnel 
and in ground effect. Thus, at low wake angles, the floor of the tunnel provides the major 
contribution to the level of interference at the model. (This situation might be altered if 
the model were placed differently in the tunnel. The correspondence would be greater if 
the model were below the centerline of the tunnel; however, if the model were placed well 
above the centerline the ceiling would assume greater importance.) At higher wake angles, 
this relationship fails (for example, 6 XjXjX in fig. 13(a)) reflecting a much greater effect 
of the walls and ceiling on the interference factors. 

If a lateral wake deflection is superimposed on the vertical deflection, as in figure 14 
(xh = 60°) and in figure 15 (xh = 30°), significant changes are observed in the interference 
factors. First, the antisymmetric group of factors are no longer zero and may assume 
substantial values (for example, ®x,y,z figs* 14(b) and 15(b)). Secondly, because the 
wake strikes the wall first (case II) for Xy > 68.95° at Xh = 60° and for Xy > 40.89° 
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at Xh ~ 30° (eqs. (44) and (46)), significant changes in the trend of several of the factors 
may be observed at these points (for example, §z,z,y in figs. 14(i) and 15(i)). No simi- 
lar alteration is observed in the ground-effect interference factors since there is no wall 
present to alter the wake pattern in ground effect. In general, it will be observed that 
there is some increase in almost all of the interference factors as the horizontal skew 
angle decreases. 

Effect of XH*“ Similar presentations of the interference factors in figures 16 to 18 
directly illustrate the variation of the interference factors with xh* Figure 16 has been 
prepared for Xy = 90° where the wake never reaches the floor. Figure 17 treats a ver- 
tical skew angle of 60° where the wake strikes the wall first for all Xh < 49.11°, and fig- 
ure 18 treats xy = 30° where the wake strikes the wall first for all Xh < 21.05°. In 
ground effect, varying xh produces relatively mild changes. This would be expected 
since with suitable resolutions of the axis system, the same results could be obtained from 
the Xh = 90° calculations. There is little correlation between the interference factors 
for ground effect and for the wind tunnel since the predominant effects on the variation of 
the factors with Xh hi the tunnel are caused by the walls in this case, and these walls are 
not present in ground effect. Some noticable changes in slope with respect to xh may be 
observed as the wake changes between striking the wall first and striking the floor first 
(for example, 5 x ,x,x in fig* 18(a)). Note that in figure 16 where xh ~ 90°, the antisym- 
metric group of factors with respect to z (those in which the subscripts contain z an 
odd number of times) are zero in the wind tunnel. They are not zero in ground effect 
because the boundaries are not symmetric with respect to the X-Y plane. 


Distribution of Interference Factors Over the Principal Axes 

Longitudinal axis .- Figures 19 to 24 show the distribution over the X-axis of the fac- 
tors representing the interference velocities for vertical skew angles of 90°, 60°, and 30° 
in combination with horizontal skew angles of 90° and 60°. 

The physical concepts of nonuniform interference as rotational rates require that the 
interference velocities be simple linear functions of distance along a given direction. If 
this relationship is valid, at least approximately, in the region occupied by the model, it 
may be possible to obtain relatively simple closed-form expressions for the corrections 
resulting from the nonuniformity of the wall effects. Figures 19 to 24 examine the degree 
of approximation involved in such a linearization by comparing the values obtained by a 
direct calculation at a series of points along the X-axis with those obtained by using the 
interference gradients in the linear relationship 


5 i,j ~ 6 i,j 


r° H 



(113) 
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Except for the combinations of skew angles which require certain factors and their 
gradients to be zero, it will be observed that very significant differences are obtained 
between the values at the model and the values corresponding to other points along this 
axis. These differences will affect measured pitching moments, particularly if the model 
has a long tail length. In most cases, the linear relationship using the interference gradi- 
ents (eq. (113)) yields a reasonable approximation to the directly calculated interference 
factors for a distance on the order of one-third of a semiheight along the axis. This 
usable distance may shrink dramatically for wakes which have large deflections in either 
direction. A noticeable example of this trend is 6 ZjX (figs. 21(c) and 24(c)), where the 
model must be relatively small longitudinally if equation (113) is to be used. 

Lateral axis .- A similar comparison over the lateral axis is given in figures 25 
to 30. Observe that the y-distance is nondimensionalized with respect to the semi- 
width B in these figures. Thus, the lateral analog of equation (113) becomes 


6 i,j = 


H 


+ y § 5 W 


^=0 

H 


(114) 


The interference distribution over the Y-axis is far from uniform. For the inter- 
ference factors which are antisymmetric with respect to y, the linear relationship given 
by equation (114) will yield reasonable results over about the central quarter of the tunnel. 
The situation is far worse with respect to the symmetric interferences. The lateral varia- 
tion of these interferences may be totally lost if represented by equation (114). (See, for 
example, 6 ZjX and 5 ZjZ in fig. 27(c).) These interferences, being symmetric, cannot 
be represented as equivalent rates of rotation. Instead, their interaction with the model 
must be studied in detail. The effect on the measured data may be only a change in stall 
angle for an unswept wing; however, these symmetric interferences can produce signifi- 
cant pitching moments on a swept wing. 

Vertical axis .- A similar comparison along the vertical axis is presented in fig- 
ures 31 to 36. In this presentation, the interference factors are obtained from the gradi- 
ents as 


6 i,j _ 5 i>j 


2U 0 + I 6 ‘ d ’ Z 

H 


^=0 

H 


(115) 


As might be expected, the behavior over this axis is similar to that over the Y-axis 
since both axes are transverse to the tunnel stream axis. The anstisymmetric factors 
with respect to z are represented reasonably well by the linear relationship for about 
the central quarter of the tunnel height; however, linear representation of the symmetric 
interference components is satisfactory. The rapid increase of interference near the floor 
when the wake is sharply deflected might be expected to lead to a reversed flow at the floor 
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itself as in references 8 and 11 with limitations on the maximum wake deflection (or mini- 
mum forward speed) which can be tolerated in the wind tunnel (refs. 2, 10, 13, and 14). 


Effect of Finite Span 

The preceding results, and the theoretical treatment given in the present paper, are 
all obtained by assuming that the model is vanishingly small. This artifice considerably 
simplifies the mathematical treatment and yields reasonably accurate interference factors 
for models which are quite small with respect to the wind tunnel. In practice, models tend 
to be rather too large for such a simplistic treatment, and the interference factors may be 
significantly altered for models of different sizes or of different configurations (refs. 1, 3, 
5, and 6). In such cases, the theoretical treatment of the vanishingly small model can be 
used as a building block to obtain the appropriate results (by superposition) for models of 
large size. Reference 6 presents a systematic means of accomplishing these superposi- 
tions, and reference 7 provides specific FORTRAN programs (based on the theory of 
ref. 5) for a wide variety of configurations. The identical treatment has been applied 
herein to develop interference factors for an arbitrary swept wing of finite size. Since the 
one significant difference between the present programs and those given in appendixes B, 

C, and D of reference 7 is the provision of an arbitrary yaw angle, these programs are not 
reproduced herein; however, a few sample calculated results are now presented. These 
results should indicate the nature and magnitude of some of the effects of finite size. 

Effect at model .- Figures 37 to 40 illustrate the effect of the ratio of model span to 
tunnel width on the interference factors. In all cases, the model is assumed to be an 
unswept wing, centered in the wind tunnel, and mounted at zero angle of yaw. While there 
are a few exceptions (for example, in fig. 38(g) for 6 ZjX ,x when Xy > 35°), as a general 
rule the interference factors decrease continuously as the span-width ratio a increases. 
This should not be construed as a decrease in the interference, for the actual interference 
depends on the product of the interference factor and the area ratio 


Am _ 7rs2 _ t r /s\2b __ n 2 V 
At 4BH 4 \B/ H 4 7 


(116) 


where A m has been taken as the area of a circle circumscribing the wingtips. Exami- 
nation of figures 38 to 49 indicates that the reduction in the interference factors is seldom 

A yy. 

of sufficient size to overcome the large increase in the area ratio ^r» Thus, the inter- 
ference in the tunnel will increase with increases in model size, but not quite in proportion 
to <j2. The effect of a is sufficiently great that finite span should certainly be included 
in the calculations whenever the model span exceeds about 10 percent of the wind-tunnel 
width. 
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Effect on distribution over the principa l axes.- The distributions of interference fac- 
tors over the principal axes have been obtained for a centered unswept wing with zero yaw 
when the wing spans half of the tunnel width (a = 0.5). These distributions are presented 
in figures 41 to 58. The directly calculated distributions are compared with those obtained 
from the gradient using the values calculated at the center of the wing and also using the 
values obtained by averaging over the span of the entire wing. The lateral distributions 
(figs. 47 to 52) are presented in terms of y/s where s is the semispan of the wing. 

In terms of this parameter, equation (114) may be rewritten as 


6 i,j " 6 i,j 


y^ 0 +aY s 5i ’hy 

H 


H 


(117) 


Examination of figures 41 to 58 indicates that the extent of agreement between the 
direct and linearized 6^ j is about the same as for the vanishingly small model. Along 
the lateral axis there is actually some small improvement presumably because the wake 
itself is spread over this axis. As might be expected, using the average values of the 
slope and intercept in equation (117) yields a better approximation to the distribution over 
the Y-axis; however, there is still a substantial loss of information if the symmetric fac- 
tors (with respect to y) are represented by the linearized form. The values along the 
X- and Z-axes are represented better by using the centerline values in equations (113) 
and (115) since the directly computed values are not averages over a finite span but repre- 
sent only the value on the axis itself. 

Effect of differences in configuration .- The effect of differences in configuration 
is briefly examined in figures 59 to 62 for two combinations of wake skew angles 
( X H ~ Xy = 90 ° and Xf£ = Xy = 60°). As before, the span of the wings is half the width 
of the tunnel, and the apex of the wing lifting line is centered in the tunnel. Only the fac- 
tors yielding the velocities and their gradients in the X-direction are presented. The con- 
figurations compared have either 0° or 45° of sweep and a yaw angle of either 0° or 45°. 

In all cases, the distribution shown is along the lifting line of the wing. 

These figures indicate substantial differences in the gradient interference over the 
span. Conversion of these interference velocities into their effect on the span loading 
would result in corrections to both pitching and rolling moments. Observe that the gradi- 
ents in the X-direction are often thought of as an effective camber, but that for the 45° 
swept wing with 45° of yaw the effect on one half of the wing is more nearly of the nature 
of an altered dihedral angle since that side of the wing lies along the X-axis. 

Even though the effects of interference distributions such as these can be large, it is 
somewhat difficult to envision complete corrections being made in a routine fashion. Even 
though reasonably complete treatments of spanwise loading are available for wings, it 
would appear that the required expenditure of effort and computer time would be excessive. 
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It might be considerably more effective to examine the interferences prior to test and to 
choose a model size such as to minimize the corrections to the point where cruder approx- 
imations to the actual corrections may suffice. Such a pre- examination will most likely 
indicate that models intended for lateral-directional stability testing should be somewhat 
smaller than models intended solely for symmetrical tests (such as performance or longi- 
tudinal stability). 


CONCLUDING REMARKS 

A theory has been presented which predicts the interference velocities caused by the 
walls of a wind tunnel as well as the gradients of these interference velocities over the 
principal axes of the tunnel. The theory allows for large wake deflections in both the lat- 
eral and vertical directions. It includes available V/STOL interference theory, where the 
wake is deflected only in the vertical direction, and conventional interference theory, where 
the wake is completely undeflected, as special cases. Various symmetry and interchange 
relationships are developed which by themselves, provide significant restraints on the 
interferences, and a large number of numerical results typifying the behavior of the inter- 
ference factors have been presented. The equivalent results for ground effect appear as 
a degenerate case of wall interference. 

Although the application of this interference study to the correction of data has been 
discussed in general terms, its immediate application depends upon the availability of ade- 
quate theoretical aerodynamic treatments of the effects of a nonuniform flow field on the 
model. Even for conventional configurations, where numerical vortex- lattice theory could, 
in principle, be used, such corrections may possibly be uneconomically time consuming. 
For more unconventional aircraft, adequate theoretical treatments often do not exist and 
rigorous corrections are not possible. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 6, 1971. 
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APPENDIX A 


PARTIAL DERIVATIVES OF £ 


Let 


and 


R 0 = 


C HV 



= yjl - COS' 


'X H COS Xy 


General Case 

A 

From equation (29b) and the above definitions, -g- is 


£ = r o c hv ~f|) sin % sin X V + f h) sin X V cos X H + 
From equation (Al.l), the partial derivatives of the first order are 



sin Xjj cos x v 

(Al.l) 



- sin Xjj sin Xy 




+ sin Xy cos x H 


(A1.2) 


(A1.3) 




+ sin x H cos Xy 


(A1.4) 


From equations (A1.2) to (A1.4), the partial derivatives of the . second order are 


° 2 (r) _ £hvJ 

[ (^S) 2 ] 

0/x\ 2 R 0 

L 

w 


d2 (w) C HV 

[Jtfl 

Ro 

«o 2 J 


(A1.5) 


(A1.6) 
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APPENDIX A - Continued 
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W “° 
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(A1.7) 

(A1.8) 

(A1.9) 

(A1.10) 


From equations (A1.5) to (A1.10), the partial derivatives of the third order are 



- 1 



- 1 


(Al.ll) 


(A1.12) 


(A1.13) 


(A1.14) 


(A1.15) 


(A1.16) 
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APPENDIX A - Continued 



Special Case of Xy = 90° 

Setting Xy = 90° in equations (Al.l) to (A1.20) yields 

£ = R o - £ g) sin X H + (? |) cos % 

The partial derivatives of the first order are 



The partial derivatives of the second order are 



(A1.17) 


(A1.18) 


(Al.l 9) 


(A 1.20) 


(A2.1) 


(A2.2) 


(A2.3) 


(A2.4) 


(A2.5) 
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APPENDIX A - Continued 



(A2.6) 


(A2.7) 


(A2.8) 


(A2.9) 


(A2.10) 


(A2.ll) 


(A2.12) 


(A2.13) 


(A2.14) 


(A2.15) 
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APPENDIX A - Continued 



Special Case of Xjj = 90° 

Setting x H = 90° in equations (Al.l) to (A1.20) yields 



(A2.16) 


(A2.17) 


(A2.18) 


(A2.19) 


(A2.20) 


(A3. 1) 


(A3. 2) 


(A3. 3) 


(A3. 4) 
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APPENDIX A - Concluded 


The partial derivatives of the second and third orders are as given in equations (A2.5) to 
(A2.20). 


Special Case of Xjj = Xy = 90° 
Setting Xji = Xy = 90° in equation (Al.l) yields 


A 

h 


= R o- 



(A4.1) 


The partial derivatives of the first order are 



(A4.2) 


(A4.3) 


(A4.4) 


The partial derivatives of the second and third order are as given in equations (A2. 5) to 
(A2.20). 
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APPENDIX B 


PARTIAL DERIVATIVES OF 

h 


In this appendix, the equation numbers are chosen to correspond with appendix A. 


General Case 


A» 

From equation (40), — is 


so that 


so that 


- Rq^HV " s ^ n s * n X V - h) S * n X V C0S X H + (j* h) S * n X H COS X V 


8 (y) h)°hv 

W R ° 


- sin Xy COS X 


H 


(Bl.l) 


(B1.2) 


Special Case of Xy = 90° 
TT = R ° " (f ll) sm X H “ (f h) cos X H 

»MM. cosx 

W R ° H 


(B2.1) 


(B2.2) 


All other values of the partial derivatives of -j— are identical with those in 


appendix A. 
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APPENDIX C 


INTERFERENCES AND INTERFERENCE DERIVATIVES AT THE 
CENTER OF THE LIFTING SYSTEM IN GROUND EFFECT 


Notation 

Because of the length of many of the expressions derived herein, the following abbre- 
viations are used both in this appendix and in table I: 


c H = cos x H sh = sin x H 

cy = cos x v s y = sin x v 

Cjjy = \/l “ COS Xjj COS Xy 


Method of Derivation 

The interference factors for the velocities are obtained by a straightforward substi- 
tution of and j into equation (65). Similarly, the interference factors for the 

derivatives of the velocities are obtained from equation (67). The values of A and A' 
and the derivatives thereof are given in table I. Substantial manipulation of the trigono- 
metric terms is required to obtain the final results. 

The symmetries discussed earlier (eqs. (70), (71), and (72)) are of substantial value. 
One velocity interference factor and 12 interference factors for the gradients are obtained 
directly from the remaining factors. 

The expressions presented herein are valid only for Xy = 90°. Simpler expres- 
sions for Xy > 90° may be derived from equations (68) and (69) in the main text. 


Interference Factors for Velocities Caused 
by Forces in the X- Direction 
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APPENDIX C - Continued 


When Xfj is 90 c 


5. 


-I 4s c 2 - 3s 2 c 2 , 1 / V ^ 
:,x - if r s v c v Js v c v + 2 i + c v J 


(Clb) 


Equation (Clb) is identical, except for notation, with the equivalent expression given 
as equation (A16) of reference 5. 

When Xy is 90° 

. 2 
"H 

w 


6 


x,x 


When x H and Xy are 99° 


5x,x = 0 


(Clc) 


(Cld) 


x,y*' 
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c hv 3 ( c 
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s 2 C 2 2s 4 c 4 
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2 C HV " S H C V 


(C2a) 


When Xjj is 90° 


6 *,y “ 0 


When Xy is 90° 

6 . 


S H C H 


x >y 2n 


(C2b) 


(C2c) 


J x,z-' 


X,y 7T 



(C3a) 


When x H is 90° 


6 x,y = 7 r ^ s V c y 3 " 4c V 3 " S V C V ~\ taxi 2 


(C3b) 


which is identical, except for notation, with the equivalent expression given as equa- 
tion (A12) of reference 5. 
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APPENDIX C - Continued 


When x v is 90 c 
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When both Xjj and Xy are 90° 
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Interference Factors for Velocities Caused by Forces 
in the Y-Direction 

6y,x- See equations (C2). 
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When both Xjj and Xy are 90° 
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(C3d) 


(C4a) 


(C4b) 


(C4c) 


(C4d) 


(C5a) 
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APPENDIX C - Continued 


When Xjj is 


90° 


6 y,z - 0 


When Xy is 90° 




When Xjj is 


7 ZTI 

Interference Factors for Velocities Caused by 
Forces in the Z-Direction 
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Equation (C6b) is identical, except for notation, with the equivalent expression given as 
equation (A8) of reference 5. 

When x v is 90° 
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When both Xjj and Xy are 90 c 
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APPENDIX C - Continued 


When x v is 90° 
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When x H is 90° 
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Equation (C8b) is identical, except for notation, with the equivalent expression given as 
equation (A4) of reference 5. 

When x v is 90° 


6 z> z " 2 u 


When both x h and Xy are 90° 
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Interference Factors for the Derivatives of the Velocities 
Caused by Forces in the X-Direction 
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s C 4 
Sy^HV 

3s h 3c v 3 ( c hv 2 " s h 2s v 2 ) 

> 

a 

to 

1 

_ 1 

s v 3c hv 2 



(C8c) 


(C8d) 


(C9a) 
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APPENDIX C - Continued 


When x H is 90° 
6 x,x,x = 

When Xy is 90° 
6 x,x,x = 
When both Xjj and 
^x,x,x = 
5 x,x,y 

6 x,x,y = ' 


When x H is 90° 

5 x,x,y = 0 
When Xy is 90° 

6 x,x,y = - 


10Sy c y5 - 7 s T7 -c T 7 .3 _ 12c Tr ^ + 8c Tr 3 - -1 


+ 2c 


V V 


V 4 1 + c- 


^)tan$ 


(C9b) 



Xy are 90° 

1 

~ 3i r 


(C9c) 


(C9d) 


C H 

2g 3 g c 3(r 2 _ c_ 2„ 2^\ 

S V C V \THV 0S H s V ) 

4s 3 c 3 (c 2 - 3s 2 s 
S H C V THV JS H S V ) 

7T 

c ® 

^HV 

C 5 
^HV 


8S H 5c y 3 SyC H y ^ S R 2 Sy3 

C HV 3 4 ( C HV _s H c v) 2 2 ^HV" s H c v) 3 

s h 3c v 3 ( c hv 2 ~ 3s h 2s v 2 ) [ s h 3c v 3 ( c hv 2 " s H 2s v f 

s 3 c 2 

V HV b y 

2s 2 3 

^ S H C V ^ S H C V 

S V C HV 2 S V 3 


(ClOa) 


(ClOb) 



(ClOc) 
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APPENDIX C - Continued 


J x,x,z*‘ 


1 

,A,Z, 


^x.x.z - 


2 s H 4 c V 4 / 2 „ 2 c 2 \ 4 s h" c v’J 2 s H% 1 ^ C HV 
r ( 5S H S V - C HV ) ~ c B {Cmr + s 


"HV 


4 " 4 ( 3 s H 2 s y 2 " C HV 2 ) _ 5 f] 


"HV 


( C HV + S v) 

4 s h 4 c v 4 ( c hv 2 " s h 2 s v 2 ) + s h 4 c v 4 ( c hv 2 " 3 s h 2 s v 2 ) 
C HV 4 ( C HV + s v ) 2 s v 2 c hv 4 


c 2p _ q 3 r 2s 6 C 4 
C H ^HV S H C V H C V 


2 ( C HV " s H c v) s v 2 c h V S 


When x H is 90° 


6 X,X,Z ~ (lOSy® “ 12 Sy^ - 19 Sy 4 + 24 Sy 3 + 9 Sy 2 - 12 Sy 1 + 2 C y + 1 


When Xy is 90 


c 2 

5 =fS_ 

3 x,x,z 2 


When both x H and Xy are 


90° 


^X,X,Z - 0 


®x,y,x ,_ See equations (CIO). 


J x,y,y* 


5 =^a 

x,y,y 7 t 


2 s h 3 s v c v 3 


(s R 2 - 4 s y 2c H 2 ) 4 s h 3 c v 3 ( s h 2 ~ 2 s V 2 c H 2 ) 


C 6 
l 'HV 


HV 

s h 3 c v 3 ( s h 2 ‘ 2 s v 2 c h 2 ) . 8 s h 3 c h 2 c v 3 2 s h 3 c h 2 c v 3 

' A O " O 


S V C HV 4 


"HV 


s C 2 
S V^HV 


G 5p 3 3p 2p 3 s C 

S H C V 4S H C H C V S V L 'HV 

+ « ' \2 + 


s 3 C 2 
S V C H 


o 3 C 2 

s y ^HV 


’V 


4 ^HV “ S H C v ) 2 2 (pHV " S H C v)" 


. 4 
HV 


(Cl la) 


(Cllb) 


(Cllc) 


(Cl Id) 


(Cl 2a) 
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APPENDIX C - Continued 


When x H is 90° 

6 — 1 

°x,y,y - i 

When Xy is 90° 

6 x,y,y = " i f (* ~ 2 c h 2 ) 

When both x H and Xy are 90° 


|S V C V 3 - s V 3c V - 4c V 3 - 4^ + c^ 


tan 


Xy 


g L 

°x,y,y “4^ 


6 


%y,y m 


°x,y,z 


S H C H 


1 Qc 4 P 4 q 2 fie 4p 4e 3s 4 r 4 Oe 4p 4 

1U&JJ Ly by Objj Cy by fibjj Ly ZSjj Cy 

- — + “ — + 


C 6 
'"'HV 


C 5 
U HV 


C 4 cs 2c 2 
’■'HV S V '-'HV 

8s 4 

H "V "V os H C V 

C HV 5 ( C HV + s v) C HV 4 ( C HV + S v) 2 C HV 3 ( C HV + s v) 

“I 


4s ^s ^ 
’®H “V C V 


4s 4 2 4 

^°H S V C V 


s 2 
S V 


When Xfj is 90° 
5 x,; 

When Xy is 90 1 


2 ( C HV ‘ s H c v) 2 


6 x,y,z - 0 


* _ S H C H 

x,y,z " 277 

5x,z,x* - See equations (Cll). 
6 x,z,y*“ See equations (Cl 3). 


6. 


X,Z,Z‘~ 


6 = fH 

%Z,Z 77 


2s h 3s v c v 3 ( s v 2 " 4 s h 2c v 2 ) 4s h 3c v 3 ( s v 2 ~ 2s h 2c v 

C 6 C 5 

'-'HV '“'HV 


s h 3s v c v 3 2s h 5c v 5 s h 3c v 3 1 


C ^ 
'■'HV 


2s H c V t S H C V | 1 S V 
S V C HV 4 + S V C HV 2 + 2 CHV " SRCV _ 


(Cl 2b) 

(C12c) 

(C12d) 

(C13a) 

(C13b) 

(C13c) 

(C14a) 
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APPENDIX C — Continued 


When x H is 90° 


6 x,z s z = i ( 12 c V 5 - 10 s V c V 5 - 4 c V 3 + 5 s V c V 3 + S V C V + \ 1311 ^T) 
When x v is 90° 


* _ S H 

°x,z,z - 27 


When both Xjj and Xy are 

* 1 

°x,z,z - ^ 


90° 


Interference Factors for the Derivatives of the Velocities 
Caused by Forces in the Y- Direction 
5y,x,x »~ See equations (CIO). 

6 y,x,y“ See equations (CIO). 

See equations (C13). 

5 y,y,x*" See equations (C12). 


6 


(C14b) 


(C14c) 


(C14d) 


y>y>y-' 



^ S H^ S V C V^ 

( 2s y 2c H 2 

- 3s h 2 ) 

12 s h^ c V^ 

Qo 3 p 2 C 3 
0S H C H C V 

3s H 

5 C 3 

c v 

TT 


C HV 6 


+ c 5 

^ HV 

C 3 
^HV 

s C ’ 
Sy'-HV 

Qq 5 r 3 

6s h c v , 

4s h 3c h 2c 

v! + 3. 

Sy C HV 

c 3 p 2 

S V C H 


(Cl 5a) 

s v 3c hv 2 

s 3 
s y 

5 1 

( C HV " s H c v) 2 

2 ( C HV " s H c v) 3 


When x H 
6 . 


is 90° 
= 0 


y,y,y 


(Cl 5b) 


When 


Xy 


is 90° 


5 --Sfe 2+i 

°y,y,y ~ 2n f h + 2 


(C15c) 
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APPENDIX C - Continued 


y»y>z-~ 


"y 5 y,z - tt 


ir 2 s H 4 c V 4 ( 4 s V 2 c H 2 - S H 2 ) 8 s H 4 s V c H 2 c V 4 4 s H 6 c V 4 


c hv 5 ( c hv + s v) 


1 4 s H 6 c y 4 8 s h 4 c h 2 c v 4 | 2 s h 4 c h 2 c v 4 s h 6 c v 4 

c hv 4 ( c hv + s v ) 2 c hv 3 ( c hv + s v) c hv 4 s v 2 c hv 4 


H "V 


HV I HV T °Vj HV rHV T D v 


’V HV 


, 2 s H 4 c H 2 c V 4 t C HV s V 2 c H 2 

s v 2 c hv 2 2 ( c hv s h c v) 2 ( c hv - s H c v) 2 


’V H 


When x H is 90 c 


5 y,y,z tt f 8 " 12s v + 3 s y 2 + 4 s v 3 " 2 v 4 " \ 1 + c . 


When Xy is 90° 


y,y,z - 27 r 


When both x H and Xy are 90° 


y,y,z 277 


^y,z,x* - See equations (C13). 


5 y>z,y" See equations (C16). 


J y,z,z- 


c H 2s H 3s v c v 3(4s H 2c v 2 - s^j 4s H 3c v 3 (2s H 2 c v 2 - S y 2 

y>z,z 77 p 6 ■ p 5 

U HV C HV 


s h 3 s v c v 3 


1 s v 


9 c 5p 5 o 3r» 3 
ZS H C V S H c v 


^ ^HV S H C V s y^HV^ S V^HV^ 


ii i in ll l mil him i mm 1 1 i i i i II I ■ 



APPENDIX C — Continued 


When Xjj is 90° 


6 y,z,z = 0 

(C17b) 

When Xy is 90° 

6 - _ fa 

°y>z,z - 2 tt 

(C17c) 


Interference Factors for the Derivatives of the Velocities 
Caused by Forces in the Z- Direction 


5, 


z,x,x*‘ 


6, 


z,x,x - ^ 


i° s h 6s v 2c v 4 s h 4c v 4 ( 5 " 3c h 2 ) s h 4c v 4 ( 1 " 2c H 2 ) 


C 6 
'-'HV 


C 4 
'-'HV 


s 2c 2 
S V '-'HV 


'HV 


a 2 S 2 

S H S V 


i(CHV-s H c V ) 2(C HV -S H C V )2 


(Cl 8a) 


When Xjj is 90° 


5, xx =; (lOc,.. 6 _ 5c V^ “ C V^ “ ^ i 

Z,X,X 7T l V V V 21 + Cy 


When V-Tr is 90° 


6 


X v 

C H 2 
Z,X,X 27 T 


(Cl 8b) 


(C18c) 


When both Xjj and Xy are 90° 
5z,x,x = 0 


(C18d) 


6 


z 5 x,y 


6 


S H C H P S H 4s V 2c V 4 3 s H 4c V 4 


z,x,y 


C 6 
'-'HV 


2s h 4c y! * I c av 1 s h c v) (cl9a) 


'HV 


Sy 2c HV 2 2 C HV ■ S H C V/ 
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APPENDIX C - Continued 


When Xjj is 90° 

6 z,x,y = 9 (C19b) 

When Xy is 90° 

6 z,x,y 2tT~ (C19c) 


J z,x ? z*' 


6 = £5 

U Z,X,Z 77 


2s H 3s v c y 3 ( s v 2 ~ ^ S H^ C V^) 2s H^ c V^ 2 s H^ c V 3 


'HV 


S V C HV 


s C 2 ‘ 
s y^HV 


_ s h 3s v c v 3 _ 1 s v 

C H y 4 " 1 C HV - S H C V 

When Xjj is 90° 

^Z,X,Z = “ ft ^-9SyCy9 + SyCy 2 + SyCy + ^ tan ~j 

When Xy is 90° 

5 =-M 

Z,X,Z 277 

When both Xjj and Xy are 90° 

6 - - — 
z,x,z 2 t 7 


6 z,y,x* _ See equations (C19). 

6 z,y,y" 


(C20a) 


(C20b) 


(C20c) 


(C20d) 


5, 


z,y,y " 77 


2s h 4c v 4 ( 4s v 2c h 2 ~ s h 2 ) _ 2s h 4c h 2c v 4 s h 6c v 4 _ 2s h 4c h 2c v 4 


C 9 
^HV 


'HV 


s v 2c hv j 


c 2c 2 
s y '-'jjV 


'HV 


s 


2 C 2 
V C H 


2 ( C HV “ s H c v) 2 ( c H y " s H c v) 2 


(C21a) 
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APPENDIX C - Continued 


When x H is 90° 

6 z,y,y = - ?^ C V 4 + C V 2 + 1 ■ \ 1 + C v) 
When Xy is 90° 

6 =. fs! 

°z,y,y 2ir 

When both Xjj and Xy are 90° 

« _L 

°z,y,y - " 2u 


(C21b) 


(C21c) 


(C21d) 


5 z,y,z-" 


6 = ^h 

z,y,z n 


2s c ^ 

®V v 


( 4s H 


2 2 2 
c y - s y 


C ^ 
^HV 


s c ^ 
S H S V C V 

C ^ 
U HV 


2s 5 C 5 


H_V 


Stt-C 


V HV 


+ s h 3c v 3 + 1 s v 

SyC H y 2 2 C HV - S H C V 


(C22a) 


When Xjj is 90° 

\y,z - 0 

When Xy is 90° 

6 

z,y,z 2 it 

6 z,z,x*“ See equations (C20). 
6 Z)Z ,y - See equations (C22). 


(C22b) 


(C22c) 


77 
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When x H 


When x v 


When both 


APPENDIX C - Concluded 


/ /s 6p 6 30 2p 4 

6 = I f H C V ^ S H S V C V 1 

z > z > z in p 6 ~ p . 6 2 

U HV U HV 


is 90° 


6 z,z,z - J (* 0c v 6 " 3c y 4 + 2 


is 90° 

6 = — 

Z ? Z ? Z 27T 

X H and x v are 90° 

5 = _L 

z > z > z 27T 


(C23a) 


(C23b) 


(C23c) 


(C23d) 
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APPENDIX D 


FORTRAN PROGRAM TO CALCULATE THE INTERFERENCE FACTORS FOR STABILITY 


TF STS OF A VANISHINGLY SMALL MODEL IN A CLOSED RECTANGULAR TEST SECTION 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN. VERSION 2.1. TO RUN CN CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BF REOUI RED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BFFN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE FQUIVALFNT OF APPROXIMATELY 15 DECIMAL OIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY RFOUIRE MOD IF EC AT ION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF FQUAl ACCURACY. 


INPUT IS REQUIRED AT LINE ID77) AND EITHER LINE ID 891 OR ID 92). THE 
FIRST CARD MUST STATE. IN NAMELIST FORMAT. THE HORIZONTAL ICHIH) AND VERTICAL 
ICHIV) SKFW ANGLES FOR WHICH THE INTERFERENCE FACTORS ARE REQUIRED. AS MANY AS 
11 HORIZONTAL AND 10 VERTICAL SKEW ANGLES MAY BE SUPPLIED. AND THE OUTPUT WILL 
RF FOR All VFRTtCAL SKEW ANGLES AT EACH HORIZONTAL SKEW ANGLE. THE SECCNC CARD 
SPFCIFIFS UN AN A FORMAT) THE DESIRED BOUNDARY CONDITION AND THE DESIRED 
COMBINATION OF INTERFERENCE FACTORS, THAT IS: STARTING IN COl UMN 1, THIS CARD 
WILL RF AD FITHER Cl OSED TUNNEL CASE, OR GROUND EFFECT CASE: AND STARTING IN 
COLUMN 25. THF CARD WILL READ EITHER VELOCITIES ONLY. VELOCITIES AND LONGITUDIN- 
AL DFRIVATIVFS ONLY. OR VELOCITIES AND DERIVATIVES. THE THIRD CARD SPECIFIES 
IXTRA. 7 FT A . ETA, GAMMA. XOVERH, YOVERH, AND ZOVERH IN FORMAT NO. 80 (LINE 
ID 25111. IXTRA IS DEFINED IN A SUBSEQUENT LIST OF STATUS INDICATORS. THE 
REMAINING QUANTITIES ARE MERELY THE SPELLED OUT VERSIONS OF QUANTITIES DEFINED 
IN THF MAIN TEXT. AS MANY SUCH CARDS AS DESIRED MAY BE INSERTED AT THIS POINT 
PROVIDED THAT THF PRFVIOUSLY SUPPLIEO QUANTITIES ARE UNALTERED. A SET OF SAMPLE 
DATA CARDS FOLLOWS THE PROGRAM LISTING. 



PRflGR AM 

STABIL 

C I NPUT. OUTPUT. TAPE5= INPUT, TAPE6=0UTPUT I 

CD 1) 

c 




CD 2) 

c 




CD 3) 

c 

THIS 

PROGRAM 

COMPUTES INTERFERENCE FACTORS FOR STABILITY WORK 

(D 4) 

c 


AT 

A POINT NEAR A VANISHINGLY SMALL MODEL 

CD 5) 

c 




(D 6) 

c 




(D 7) 

c 


♦♦♦DEFINITIONS OF STATUS I ND IC ATQR S*** 

(D 8) 

c 




(D 9) 

c 

♦ VAR l ABI F* 

♦ ST ATUS* 

♦DEFINITION^ 

(D 10) 

c 




(D 11) 

c 

I HV90 

1 

NEITHER CHI ( H) NOR CHICV) SET TO 90 DEG BEFORE 

CD 12) 

c 



SUBST ITUTION 

CD 13) 

c 


? 

EITHER CHICHI CR CHICV) SET TO 90 DEG BEFORE SUB- 

(D 1 A ) 

c 



STITUTION 

CD 15) 

c 




( D 16) 

c 

If.HTH 

1 

CHICHI NOT SET TO 90 DEG BEFORE SUBSTITUTION 

(D 17) 

c 


? 

CHICHI SET TG 90 DEG BEFORE SUBSTITUTION 

(D 18) 

c 




CD 19) 

c 

ICHIV 

1 

CHICV) NOT SET TO 90 DEG BEFORE SUBSTITUTION 

CD 20) 

c 


? 

CHICV) SET TO 90 DEG BEFORE SUBSTITUTION 

(D 21) 

c 




( 0 22) 
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APPENDIX D - Continued 


c 

AHPR IM 

1.0 

PHI TO BE EVALUATED 

c 


-1.0 

PHI PRIME TO BE EVALUATED 

l. 

c 

I H90 

l 

CHECH) NOT EQUAL TO 90 DEG 

c 

c 

c 


2 

CH 1 1 HI EQUAL TO 90 DEG 

I V90 

l 

CHIIV) NOT EQUAL TO 90 DEG 

c 

c 

c 


2 

CHI ( V ) EQUAL TO 90 DEG 

I BOTH 

1 

EITHER CHICHI OR CHIIV) NOT EQUAL TO 90 DEG 

c 

c 

c 


2 

ROTH CHICHI AND CHICV) EQUAL TO 90 DEG 

ICASF 

1 

CASE I EQUATIONS TO BE USED 

c 

r 


2 

CASF II EQUATIONS TO BE USED 

l 

c 

I TYPE 

t 

CLOSED TUNNEL 

c 

c 

c 


2 

GROUND EFFECT 

IXTRA 

1 

ADO EXTRA CHICV) DIVIDING CASES I AND II 

c 


2 

OMITS EXTRA CHICV) 

r 

c 

c 

c 


3 

ALTER ANGLES, BOUNDARY CONDITIONS. OR REQUIRED 
INTERFERENCE FACTORS 

TCH 

1 

CHICH) NOT GREATER THAN 90 DEG OR ITYPE EQUALS 2 

c 

c 

c 


2 

CHICH) GREATER THAN 90 DEG AND ITYPE EQUALS l 

rev 

1 

CHICV) NOT GREATER THAN 90 DEG 

c 


2 

CHICV) GREATER THAN 90 DEG 


COMMON RSU80,A0VERH,CHV(2) , S I NXHC 2 ) . COSXH ( 2 I • 

1 SIN X VC 2>,COSXVC2) ,IHV90, ICHIH, ICHIV. 

2 CHORD! 3 > .PHEI2) . 1,J,K,AHPRIM 

COMMON XDVFRH.YOVERHtZOVERHtZETA, ETA. GAMMA, TAN XH. 

1 TANXV. I BOTH* IH90, IV90* ICASE.XSETl. XSET2, 

2 XSFT3.PI, IP€ 3) , IQC 3 ) , IR 1 3 > * TERM1C 2 I * TFRM2 1 2 ) . 

3 TFRM3C 2>*TEftM4<2)tTERM5l2).TERM6(2UTERM7<2) , 

4 TERM8I?) ,TE»M9|2) »TERM10(2) . TFRM1 1 < 2 ) . TERM 12 ( 2 ) . 

5 ANSFRC 2) .DELTA (3,3 ) , DELT ADC 3 • 3* 3 ) , F ACT YP ( 6 ) , ICV 
DIMENSION CHIHIU) .CHIVC I 1 ) « VHCLO ( 1 1 ) * RUNTYP I 4 ) • 

1 PRTXY7( 3),PTXYZ( 9) 

DATA PRTXYZ/1HX, 1HY « 1 HZ/ , GRNT S 1 , F ACT ST/6HGR0UND, 6H / , 

1 PTXY7/4HX.X) , 4HX,Y) f 4HX,Z) ,4HY.X1 *4HY . Y ) • 4HY « Z ) • 

2 4H 7 . X ) « 4H7 * Y ) , 4HZ , Z ) / 

NAMF1 1ST /ANGLES/CHIH*CHIV 
CHV( 2)=SINXV< 2)=SINXH( 2) *1.0 

cnsx vi 2 i=cosxhi 21 = 0*0 
1 pi 1 i=ipi3)=ioc 1 )=i;ri 11=2 

IP! 2I = IG<2 )=I0<3)= IR< 2)=IRC3)=1 
PI=3. 14159265358979 
RAfl VSN= P I / l 80* 0 

49 IFIRST=0 

DO 48 1=1,11 

48 CHIHI I l=CPTVII) = VHOI DII)=0.0 
READ 15# ANGLE S ) 

IF I FOF , 5 ) 375,50 

50 AHPR I M= 1*0 
I X TR A=2 

I HV90=I CHIH=ICHI V= CH90= I V9 0= I CH= l C V= 1 
NOHOR 7= NQVERT=NOVS AV=0 


( 0 23) 
iO 24) 
( 0 25) 
Cl) 26) 
(D 27) 
(I) 28) 
(0 29) 
CD 30) 
(D 31) 
(D 32) 
( 0 33 ) 
(0 34) 
(D 35) 
(0 36) 
( U 37) 
( ) 38) 
t D 39 ) 
( D 40) 
(0 41) 
(0 42) 
( D 43) 
ID A4) 
( 0 45) 
(D 46) 
(0 47) 
(D 48) 

( 0 49 ) 
(I) 50) 
(051) 
(0 52) 
(D 53) 

( D 54) 

( 0 55) 
10 56) 
(0 57) 

( L) 58) 
(0 59) 
(D 60) 
( 061 ) 
(0 62) 

( 0 63) 
(D 64) 
(0 65) 
10 66 ) 
(0 67) 

( 0 68 ) 
(0 69) 

( 0 70) 
(0 71) 

( 0 72) 

( 0 73) 
(0 74) 
(0 7 5) 
CD 76) 

( D 77) 

( D 78) 

( D 79) 
(0 80 
(081) 

( 0 82) 
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APPENDIX D - Continued 


Dfl 51 1 = 1,3 

nn 5i j=i * 3 

DFL TAf I ,J)=0.0 
DO 51 K = 1 » 3 
51 DFLTAOC I,J,K)=0.0 

IF( IFIRST.EO.O) GO TO 90 

RFAD <5,801 I XTRA ,Z ETA, ETA, GAMMA* XOVERH, YOVERH, ZOVERH 
IF ( I XTRA* EG* 3 1 GO TO 49 
I FI FOF, 5 1 375,109 

90 RFAD <5,1001 R UNTYP, F ACTYP , I XTR A . ZETA, ETA , 

1 GAMMA, XOVERH, YOVERH, ZOVERH 
IF (EOF, 51 375,101 

101 IFIRST=1 
10? DO 105 M=l, 1 1 

105 VHOL DIM )=CHI V CM) 

106 I F < R UNT YPC 1 l-GRNTST I 108,107 f 108 

107 I T YPE=? 

ZFT A = l • 

FT A= 1 • 

GAMMA=?. 

GO TO 109 

108 l TYPF= 1 

109 A 7FT A=Z FT A 
AFTA^ET A 

on no m=i * io 

IFICHIHC Ml) 130,130,1 10 

110 N0H0R7=N0HQRZ+1 
130 00 140 M=1 , 10 

IF< CH IV <M) 1160,160, 140 
140 MOVER T= NO VERT ♦ l 

160 NOVS A V=NOVERT 

00 365 T ANGLH=1 , NOHORZ 
AC HI H=C HIH C I ANGLH1 
NOVFR T= NOVSAV 
IF IITYPE.EO.?) GO TO 161 

WRITF 16,9031 GAMMA, AZETA , AETA, XOVERH, YOVERH, ZOVERH 
GO TO 16? 

161 WRITF (6,904) XOVERH , YOVERH, ZOVERH 
16? I CH= 1 

IF ( CHI HI IANGLH) .LE.90..0R. ITYPE.E0.2) GO TO 164 
ICH=? 

FT A=? . O-AFT A 

C H I H ( IANGLH ) = 1 80.0-ACHIH 
YOVF RH=- YOVERH 

164 00 165 M= 1,11 

165 CH I V ( M ) =VH0L0 ( M ) 

IF ICHTHI IANGLH). EO. 90, 01 GO TO 2?0 
170 STNXHI 1 1=SINICHTHI I ANGLH 1 *RADVSN I 
CnSXHIl l=COSCCHIHC l ANGLH l*RADV SN ) 

T ANXH=S TNXH( 1 l/COSXHI 11 
l H90= 1 

IF ( ITYPF.F0.?.0R.IXTRA.EQ.2> GO TO ?30 

C 

C FIND ANO INSERT THE VERTICAL SKEW ANGIE FOR WHICH THE 

C WAKF GOES DIRECTLY TO THE CORNER OF THE WIND TUNNEL 

C ( SHARFO BY CASES I AND II) 

C 

VMTO=ATAN(TANXH*ZFTA4‘GAMMA*(2.0-ETA) l/RAOVSN 
DO 185 H= 1 • NOVERT 

IF (ABSIVMIO-CHIV(MI).LT, 0.0051 GO TO 230 


(D 83) 
(D 84) 
(0 85) 
(D 86) 
(D 87) 
CD 88) 
(D 89) 
ID 90) 
(D 91) 
<0 92) 
CD 93) 
CO 94) 
ID 95) 
(0 96) 
CD 97) 
CO 98) 
CD 99) 
CD 100) 
CD 101) 
CD 102) 
CO 103) 
CD 104) 
CD 105) 
CD 106) 
CD 1071 
CO 108) 
CD 109) 
(0 110) 
(0 111) 
CD 112) 
CD 113) 
CD 114) 
CD 115) 
ID 116) 
CD 117) 
CD 118) 
CD 119) 
CD 120) 
CD 121) 
CD 12 2) 
(0 123) 
CD 124) 
CD 125) 
CD 126) 
CD 127) 
CD 128) 
CD 129) 
CD 130) 
(D 131) 
CD 132) 
CD 133) 
CD 134) 
CD 135) 
CD 136) 
CD 137) 
CD 138) 
CD 139) 
CD 140) 
CD 141) 
(D 142) 


81 



n n n 


APPENDIX D - Continued 


1 85 CONTTNUF 

nn 190 M=l * NO VFRT 
IF ( VMin.LT.CHlVlMH GO TO 200 
190 CONT I NIJF 

CHI V I NOVFRT+l I =VMI D 
GO TO 715 

700 MOVI OC=NOVERT+l 
MOVF=MO VLOC—M 
on 710 M= l * HOVF 
CH I V ( MO VL OC ) =CHI VI MOVLOC-1 ) 

710 MnVinC=HCVLDC-l 
CHIVI MO VI OC)=VMID 
71 5 NOVFR T= NOVSA V + 1 
GO TO 730 
770 SINXHI 1 1=1.0 
COSXHI l 1 = 0.0 
TANXH=l O.OE 10 
CH VI 11=1.0 
TH90=7 

7 TO on 364 I ANGLV=1 * NO VERT 
ICV=1 

IF (CHIVI IANGLV1.E0.90.0) GO TC 740 
IF (CHI VIIANGLVI.LT. 90.0) GO TC 741 
I C V= 7 

IF IITYPE.E0.7) GO TO 241 
7 FT A= A 7 FT A/I7.0*AZET A— 1 • 0 1 
7nVFRH=-7nVFRH 

CHIVI I ANGLV 1=180. O-CHI V( IANGLV1 
741 SI NX VI l >=S INI CHIVI I ANGLV ) *RADVSN ) 

cnsxvu i=cns < chivi i anglv i*radvsni 

TANXV=SINXVI11 /COSXVI 1) 

CHV ( 1 1 = SORT I ( 1.0-f I COSXH ( 1 1**2 1 * I CDS XV I 11 **2 1 1 1 1 
I V90= 1 
I RnTH=l 

IFITANXH.GE.TANXV/I ZETA*GAMMA*I 2.0— ETA11 1G0 TO 260 
GO TO 750 
740 SINXV! 1 1 = 1 .0 

cnsxvn 1=0.0 

T ANX V= 1 0. OF 1 0 
CHV I 11=1.0 
I V90=7 
I BOT H=7 

I F ( IH90.NE.71 TBOTH= 1 
750 IC A SF = 7 

GO Tn 770 
760 f C A SF= 1 

770 XSFT1=7FTA*X0VERH 
XSFT7=X SFT 1 — T AN XV 

XSFT3= 7 F TA* I XOV ERH— GAMMA* 12*0- ETA) *T ANXH1 

CDOROTN ATF PERMUTATION LOOPS 

on 310 1=1*3 
on 310 J= 1*3 
on 310 K= 1*3 

TF I ABSCCHT VI T ANGLV 1 l.LT. 0.005) GO TO 364 
IF IK.FO.l.OR. ITYPE.E0.2) GO TO 279 
IF I FACTYPI 41 .EQ. FACT ST) GO TO 310 
TF (J.NF.l .AN0.FACTYPI6) .NE.FACTST) GO TO 310 
779 00 780 T TF RM = i * 74 


ID 143) 
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280 TFRMU ITFRM)=0*0 

GO TO ( 290*300) • ITYPE 
290 CALL WAI FFF 
GO TO *10 
*00 CAI I GRNFFF 
*10 CONTINUF 

IF(ICH.EO.l) GO TO 3000 
OO 2 092 1*1*3 

on 2092 J= 1*3 
IX=IY=2 

IF ( I.FQ.2) IX=1 
IF IJ.EQ.2) l Y= 1 

OF I T A ( I *J1=| I— !•)**( IX+IY1 l*DELTAI I * J ) 

on 2092 K= l * * 

17=2 

IF (K.FQ.2) IZ=1 

2 092 DFI. T A0( I«J«K)=( (—1* ) **l IX* I Y*I Z ) ) *DELT ADI I • J* K ) 

3000 IF I ICV.FQ.1.0R.ITYPE.EQ.2) GO TO 3001 
OO 2094 1=1.3 

DO 2094 J=1 . 3 

I X= I Y=2 

IF (I.F0.3) IX=i 
IF IJ.EG.3) I Y= 1 

DFI TA( I *J )=l 1— 1 . )**( I X* I Y ) ) *DE LT A ( I, Jl 
00 2094 K= 1*3 
17 = 2 

IF (K.FQ.3) IZ=1 

2 094 OFI TAD( I . J. K I = (<-1 . ) **! I X+ I Y+1 Z 1 ) *DELT ADI I * J*KI 
CHI VI IANGI VI = 1 80.0-CHI V( I ANGLV I 
7 OVF RH= — ZOVERH 
7FT A = A2 FTA 

*001 WRITF (6*330) ACH I H ,CHI V ( I ANGL V 1 
WRITE (6*340) PTXYZ 

WRITF ( 6.350) ( ( OEL TA I I * J ) « J* l « 3 ) • I * l • 3 ) 

IF (FACTYP(6) .NE.FACTST. AND. ITYPE. EO.l) GO TO 361 
IF (FACTYPI4). EO. FACT ST .AND. ITYPE. EO.l) GO TO 364 
00 *62 J = l * 3 

*62 WRITF C 6*370) PRTXY Z 1 J ) • ( ( OEL T AO ( I . J* K ) * K= 1* 3 ) • 1 = 1 • 3 ) 

GO TO *64 

*61 WRITF ( 6*370) PRTXYZ ( 1 ) * ( ( DELT AD ( I « 1 * K ) • K= 1 * 3 ) • I = 1 * 3 ) 

364 CONTINUE 

IF ( ICH.FO.l) GO TO 365 
CH I H ( I A NGLH ) = ACH I H 
FTA= AFT A 
YOVFRH=-YOVFRH 
*65 CONTINUE 
GO TO 50 
375 STOP 

80 FORMAT ( II * F9 . 3 . 5F 1 0 . 3 ) 

100 FORMAT 1 1 0A6/I 1 • F9» 3 • 5F10. 3 ) 

*30 FORMAT ( / / 1 4X* 9HCH 1(H) = , F6. 2* 5 X * 9HCH I( V ) = *F6.2/) 

*40 FORMAT! 1 9X . 9 I 1 H ( * A4 , 8 X ) ) 

350 FORMAT! 1 X * 1 OHOEL TA I - .- ) * 9F 1 3 . 4 I 

370 FORMAT (11H DELTA (-.-.. A1 * 1H )* F 1 i . 4* 8F 13.4 ) 

90* FORMAT (1HI// 35 X* I N T ERFERENCE FACTORS FOR STABILITY WORK IN A CLC 
ISFO TUN NFL *//43 X*AT A POINT NEAR A VANISHINGLY SMALL MODEL*/// 

23 5 X* G AMM A =*F 7 . 3 . I OX* ZET A =*F6 .3* 1 1 X*ET A =*F6.3// 

3*5X*X/H =*F7 • 3*10X*Y/H =*F6 . 3* l l X*Z /H =*F6.3/) 

904 FORMAT ( 1 H 1 /// 39X* INT ERFER ENCE FACTORS FOR STABILITY WORK IN GRCLN 
ID FFFEC T *//47X*AT A POINT NEAR A VANISHINGLY SMALL MODEL*/// 
240X*X/H =*F7. 3* 10X*Y/H =*F 6 . 3 * l OX* Z /H =*F6.3/I 
FNO 
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SURRCUTINF GR NEFF 

COMMON R.SIJBO. AOVFRH.CHVI 2 ) * S I N XH ( ? ) • CDSXHl 2) . 

1 SIN XV I 7) .COSXVI2) • IHV90f ICHIH* I CHI V, 

7 COORD ( 3 ) . PHI ( 2) • [ tJtKfAHPRIM 

COMMON xnVFRH*YOVFRH f ZO VER H * ZE T A .ET A . G AMMA • TANXH * 

1 “ TANXV. I ROTH. I H90 * I V90* IC ASE .XSET1 . XSFT2 . 

2 XSFT3.PI. IPI 3 ) « IQ( 31* IR( 3 ) * TEftMl I 2 1 * TERM? I 21 . 

* TERM3I 7 1 • TERM4 I2)?TERM5(2).TERM6I2).TERM7I 2)* 

4 TFRM8I2) . TFRM9 l 2 ) * T ERMl C I 2 ) .TERMll (21 .TERM17I7). 

6 ANSFRI?) .DELTAI3.3) tDELTADI 3. 3, 3 ) . F ACT YP C 6 ) . tCV 

IF I IV9Q.EQ.7.0R.ICV.EQ.2) GO TO 200 
COORD! 1 )=XOVERH— TANXV 
CHORD I 7 ) =YOVFRH+TANX V/T ANX-H 
C00R0I3 ) =70VERH+1 • 0 
I HV90=? 

I CHI V=7 
AHPR I M= I .0 
IF<K.NF.1>G0 TO 100 
CALL PHIIJ 
100 CALI PHI UK 

DO 110 ITER M= 1.2 
110 TFRMl (ITERM)=PHI( ITERMl 
I HV90 = 1 
I CH I V= l 

IFIK.NF.1 1 GO TO 120 
CALI PHIIJ 
120 CALL PHIIJK 

DO 125 ITERM* 1*2 

125 TFRMl I I TERM )*TERM1 ( I TERM ) -PHI I ITERMl 
AHPR T M= — 1 • 0 
ICHI V=2 
IHV90-2 

cooro (2 >=-caoRO(?) 

COORD I 3 )=— COORD (31 
IFIK.NF.IIGO TO 130 
CALI PHIIJ 
130 CAIL PHIIJK 

DO 140 I TFRM= 1*2 
140 TFRM3I I TFRM)=PHII ITERM) 

I CH I V= 1 
I H V 90= l 

I F ( K . NF . 1 ) GO TO 150 
CAII PHIIJ 
150 CAI L PHIIJK 

DO 160 I TFR M= I . 2 

160 TFRM3I I TFRMl = TFRM3< I TERM )- PH I I IT ERM ) 

200 AHPR I M=— l . 0 
I HV90= 1 
I CH I V=1 

COORD (1 >=XOVFRH 
COORD! 2 )=— YOVFRH 
COORD ( 3 ) =-ZGVER»H-7 • 0 
IFIK.NF.IIGO TO 210 
CAM PHIIJ 
210 CAI I PHIIJK 

TFRM2C 1 )=PHl I l ) 

TFRM2S2 ) = PH I 12) 

DO 215 ITFRM=U2 

215 ANSFRI f TERM )*l-7.0/P I) *1 TERMll ITERMl* 

1 < (-1.0 1**1 IP( I)*IO( Jl + IRIKl I )*ITERM2I ITERM) 

2 +TERM3 I ITERM 1 ) ) 

IFIK.NF.IIGO TO 220 

DFLTAI !*JJ=ANSER(1) 

220 DEI TAOI I.J*K) = ANSERI2 1 
RETURN 
END 
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APPENDIX D - Continued 


SUBROUTINE WALEFF 
RFAI K Oft D 

DIMFNSION TFRM(l?.2l*KORDC12.3) 

COMMON RSUBO* AOVERH*CHV( 2 ) , SI NXHC 2 ) *COSXH( 2 ) * 

1 SINXVC 2 ) *COSX V ( 2 ) , I HV90 • ICHI H* ICHI V. 

2 COORDC3I.PHI(2),l,J,K f AHPRIM 

COMMON XOVERH.YQVER*H,ZOVERH,ZETA,ETA*GAMMA*TANXH* 

1 TANXV* l BOTH* I H90, I V90 , IC ASE *XSET l* XSET 2 * 

2 XSET3*PI* IPO) * IQ ( 3 J * IR I 3 ) *TERM1 (2)*TERM2(?) . 

3 TFRM3I 21 *TFRM4C2) ,TERM5<2> ,TERM612) *TERM7< 21 . 

4 TERM8I2 )*TFRM9(2) * T ERM1 0 < 2 ) * T FRM1 1 C 2 )*TERM12 C 2 I • 

5 ANSERI?)*DFLTAC3,3) » DELT AO ( 3 * 3 * 3 ) ♦ F ACT YP C 6 ) • ICV 
FQUIVAI ENCE CTFRMI It II, TERM l( 1) ) 

OATA FACTST/6H / 

PHU n -PHI(2)=0.0 

IHV90 = 1 

ICHI H - 1 

ICHIV = 1 

AHPR I M=1 .0 

COORD <11 * XSET1 

DO 130 M = 1*7 

AM = M-4 

DO 130 N “ 1*7 

AN ^ N-4 

IF < M. FO.N.ANO.N.EO.4) GO TO 130 
COORO (?) = 7 FT AM YQVERH— 4*0*AN*GAMMA) 

COORO <31 = 7 ETA* ( ZOV ERH-4 . 0 * AN 1 
IF (K.NF.l) GO TO 100 
CALL PHIIJ 

IF <FACTYP<4) .EO.FACTST) GO TO 110 
100 CALL PHIIJK 

no no i?o iov = 1*2 

120 TERM ( 1 * I D V 1 « TERM (i,IOV) ♦ PHI (IOV) 

130 CONTINUE 

SlJMMAT I CN COOPS AND COORDINATE SET-UPS 
00 1295 M = 1*7 
AM - M - 4 

00 1 795 N = 1*7 
AN = N— 4 

YSFTl = ZETAM YOVERH-4 • 0* AM*GAMM A ) 

YSFT 2 = 7 F TAM YOVERH+GAMM AM 2.C-ETA-4.0*AMM 
YSFT3 = 7ETAM YOVERH + 2 . 0*GAMM AM 2 . 0-E T A-2 . 0* AM ) ) 

7SFTI = 7 FT A* ( ZOV ERH-4* Q*AN ) 

7 SET? - Z ETA* ( ( 70VERH-4 • 0*AN ) +GAMM A* ( 2 . O-ET A ) * ( TA NXH/ TANXV ) I 
7 SET 3 = ZFT A* ( ( 70VERH— 4 *0* AN ) — GAMMA* ( 2 *0— ET A)M TANXH/ TANXV) ) 

* INITIALIZE TFRM COUNTER * 

ITFRM = 1 

* SFT-UP COORDINATES * 

I F I I C A SE # EQ. 2 ) GO TO 140 

1 X V S FT = 2 
IXHSET = 1 

KORO ( 2*1) = XSET 2 

KORO < 2*2) = YSETl+TANXV/TANXH 

KORO < 2*3) a 7SFT1+1.0 

KORD < 3*1) = XSFT3 
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KORD 


3.21 

- 

YSET2 

(D 391) 

KORD 


3.3) 


/SET 1 + 1.0 

CD 392) 

KfiRD 


4,1) 

ss 

XSET1 

(D 393) 

KORD 


4.2) 

= 

-YSET1 

(D 394) 

KORD 


4,3) 

= 

-ZSETi-2.0 

(D 395) 

KORD 


5.1) 


XSET2 

<D 396) 

KORD 


5.2) 

= 

— YSET 1— TANXV/TANXH 

(D 397) 

KORD 


5*3) 

= 

-ZSETi-1.0 

(0 398) 

KORD 


6.1 ) 


XSET3 

(D 399) 

KORD 


6.21 

= 

-YSFT2 

(D 400) 

KORD 


6.3) 


-7SFT1-1.0 

(D 401) 

KORD 


7.1 1 

= 

XSETl 

(D 402) 

KORD 


7.2) 


YSET 3 

(D 403) 

KORD 


7.31 

= 

Z SET 1 

(D 404) 

KORD 


8.1) 

= 

XSFT2 

(D 405) 

KORD 


8.2) 

= 

YSFT3-TANXV/TANXH 

(D 406) 

KORD 


8.3) 

= 

ZSETl+1.0 

(D 407) 

KORD 


9.1) 

= 

XSFT3 

(D 408) 

KORD 


9.21 

= 

YSFT2 

(D 409) 

KORD 


9.3) 

= 

ZSET1+1.0 

(D 410) 

KORD 

UO.l) 


XSETl 

(D 411) 

KORD 

( 10,21 

= 

-YSET3 

(D 412) 

KORD 

( 10,31 

= 

-ZSETI-2.0 

(D 413) 

KORD 

(11,11 

= 

XSFT2 

(D 414) 

KORD 

( 11,21 

= 

-YSET3+TANXV/TAN XH 

(D 415) 

KORD 

(11,31 

= 

-ZSETi-1.0 

<D 416) 

KORD 

(12,11 

= 

XSFT3 

(D 417) 

KORD 

< 12,21 

= 

-YSET2 

(D 418) 

KORD 

( 12,3) 


-ZSETI-1.0 

(D 419) 

oo to ; 

?00 



(0 420) 

fXVSFT 

* 1 



(D 421) 

TXHSFT 

= 2 



(D 422) 

KORD 


2.11 

= 

XSFT3 

(D 423) 

KORD 


2.21 

= 

YSET2 

tO 424) 

KORD 


2.3) 

= 

/ SET2 

(D 425) 

KORD 


3.1 ) 

= 

XSFT2 

(D 426) 

KORD 


3,2) 


YSET2 

(D 427) 

KORD 


3.3) 

= 

7SETI+1.0 

(D 428) 

KORD 


4.1) 


XSETl 

(D 429) 

KORD 


4.2) 

= 

YSFT3 

(D 430) 

KORD 


4.3) 

= 

ZSFT1 

(D 431) 

KORD 


5,1 1 

=s 

XSET3 

(D 432) 

KORD 


5.2) 

= 

YSET2 

(D 433) 

KORD 


5.3) 

= 

7SET2 

(D 434) 

KORD 


6.1) 

= 

XSFT2 

(D 435) 

KORD 


6.2) 

3 

YSET2 

(D 436) 

KORD 


6.3) 

= 

ZSETl+l.O 

(D 437) 

KORD 


7,1) 

= 

XSETl 

(D 438) 

KORD 


7,2) 

= 

-YSFT1 

(D 439) 

KORD 


7.3) 

= 

-ZSFT1-2.0 

CD 440) 

KORD 


8,1 ) 

= 

XSFT3 

<D 441) 

KORD 


8.2) 

= 

—YSET 2 

(D 442) 

KORD 


8,3) 


-ZSET3-2.0 

tO 443) 

KORD 


9.1) 

= 

XSET2 

(D 444) 

KORD 


9.2) 

= 

—YSET 2 

(D 445) 

KORD 


9,3) 

= 

-ZSET1-1.0 

tO 446) 

KORD 

(10,11 

= 

XSFT1 

(D 447) 

KORD 

(10.2) 

= 

-YSET3 

(0 448) 

KORD 

(10.3) 

= 

-ZSETI-2.0 

(D 449) 

KORD 

(11.1) 

= 

XSFT3 

(D 450) 
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KORD < 11,7) * -YSET2 

(0 

451) 


KOR (1 (11*3) = -ZSET3-2.0 

(0 

452) 


K0R0 (12*1) = XSET2 

(0 

453) 


K0R0 (17*2) = -YSET2 

(0 

454) 


KORO (17*3) = -7SET1-1.0 

(0 

(D 

455) 

456) 


* TERMS TWO* FIVE, EIGHT, ANO ELEVEN * 

(0 

(0 

457) 

458) 

700 

ITFRM = ITFRM+l 

(D 

459) 


I F ( 1 BOTH.NE • 2 I GO TO 205 

(0 

460) 


IF ( ITFRM. EQ. 11) GO TO 1295 

(0 

461) 


ITFRM = ITERM+1 

(D 

462) 


GO TO 400 

CD 

463) 

705 

THV90 = 7 

(0 

464) 


ICHIV = IXVSET 

(0 

465) 


ICHIH = IXHSFT 

(D 

466) 


COORD (1) = KORO (IT ERM , 1 ) 

(0 

467) 


COORD (?) = KORO ( I TERM, 2 ) 

(D 

468) 


COORD (3) = KORO < I TERM, 3 ) 

(0 

469) 


IF (K.NE.l ) 00 TO 210 

(0 

470) 


CALI PHIIJ 

(0 

471) 


IF I FACTYP(4).E0*FACTST) GO TC 720 

(0 

472) 

710 

CALL PHI IJK 

(D 

473) 

770 

DO 730 TDV = 1.2 

(D 

474) 

730 

TFRM (ITERM.IOV) = T ERM( I T ERH , I DV ) +PHI ( I DV ) 

(0 

475) 


THV90 =1 

<0 

476) 


ICHIV =1 

(0 

477) 


ICHIH =1 

(0 

478) 


IF (K.NE.l ) GO TO 240 

(0 

479) 


CALL PHIIJ 

(0 

480) 


IF IFACTYP(4).E0*FACTST) GO TO 250 

(0 

481) 

740 

CALI PHTIJK 

(0 

482) 

750 

DO 760 l DV = 1*2 

(D 

483) 

760 

TFRM ( ITFRM* I DV ) = TERM ( ITERM , I DV ) - PHUIDV) 

(0 

(0 

484) 

485) 


* TERMS THREE. SIX, NINE, ANO TWELVE * 

(0 

(0 

486) 

487) 


ITFRM = ITFRM ♦ l 

(D 

488) 


IHV90 = 7 

(D 

489) 


ICHIV = 2 

(D 

490) 


ICHIH - ? 

(0 

491 ) 


IF (ICASE.F0.7) 00 TO 300 

(0 

492) 


IF (IH90-FQ.7) GO TO 450 

(D 

493) 


GO TO 305 

(0 

494) 

300 

IF ( I V90.FQ. 7 ) GO TO 450 

(0 

495) 

305 

COORD (1) = KORO I ITERM, 1) 

(D 

496) 


COORO (?) = KORO (ITERM, 2) 

(0 

497) 


COORD (3) = KORO (ITERM, 3) 

(D 

49 8) 


IF (K.NE.l) GO TO 310 

(0 

499) 


CAI L PHIIJ 

(D 

500) 


IF (FACTYPI4) .EQ.FACTST) GO TO 320 

(D 

501) 

310 

CAI L PHI IJK 

(D 

502) 

370 

no 330 mv = 1.? 

(0 

503) 

330 

TFRM ( ITERM* I DV) = TERM (ITERM.IOV) + PHI(IDV) 

(0 

504) 


IF (ICASE.F0.2I GO TO 340 

(0 

505) 


ICHIH * 1 

(D 

506) 


GO TO 350 

(D 

507 ) 

340 

ICHIV = l 

(D 

508) 

350 

IF (K.NE.l) GO TO 360 

( 0 

5C9) 


CAI 1 PHIIJ 

(0 

510) 
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APPENDIX D - Continued 


IF CFACTYP44) • EQ. FACTST ) GO TO 370 
360 CALL PHI IJK 
370 nn 380 IDV = 1*2 

380 TERM 4 I TERM* I DV) = TERM CITERM.IDVI - PHI 4IDV) 

* FXIT IF TERM NUMBER IS TWELVE * 

IF 4 ITFRM.EG.12) GO TO 1295 

* TERMS FOUR* SEVENt AND TEN * 

400 I TERM = ITERM+1 

TEC ITFRM.EQ.4.0R. ITERM.EQ.7) GO TO 405 
AHPRIM - 1.0 
GO TO 410 

405 AHPRIM = -1.0 
410 THV90 * 1 

ichiv = l 

TCHIH = l 

COORD (1) » KORO < I TERM* 1 ) 

COORO 42) = KORD 4 I TERM* 2 1 
COORO 431 = KORD 4ITERM V 3) 

IF IK.NF.l ) GO TO 420 
CAM PH 1 1 J 

IF 4FACTYPI41.E0. FACT ST ) GO TO 430 
420 CA1 L PHI I JK 
430 DO 440 IDV = 1*2 

440 TERM 4 ITERM • I 0 V 1 = TERM CITERM,IOV) ♦ PHI CIOV) 

GO TO 200 

450 IF < ITFRM.NE.12) GO TO 400 
1295 CONTINUE 

1300 IF (ICASE.E0.2) GO TO 1305 
IP0WR1 = IPC I 1+IQ4J1+IRCK ) 

IP0WR2 = IPCI) 

GO TO 1310 
1305 TPOWRl = IPCIl 

IP0WR2 a IPCII+IOC J l+IRCKl 
1310 no 131 5 IDV = 1*2 

SAnA=C-2.*7FTA**C I DV+ 1 ) ♦GAMMA/P II 
SADR=T ERM4 1 « I DV I + TERM4 2 , I DV H-TERM4 3* IDVI + 

1 C 1-1. )** I POWR 1 ) * C TERMC 4* IDV l+TERMC 5* IDV) +TERM C 6* I DVD 

SADOC 4-1. I ** IP0WR2)*4TERM47,IDV)+TERMC 8. IDV)+TERM49* IDV) ) + 

1 I 4—1.1 ♦*! 104 JUIR4K) )) *4 TERM 4 10* IDVI+TERMC 11* IDV)* 

2 TFRM4 12* TDVI > 

1315 ANSERC I DV ) =SADA*C SADB+SADEl 
IF CK.NE.l) GO TO 1320 
DFI.TA H*J ) = ANSER 41) 

1320 DFLTAD CT*J*K) = ANSER (2) 

RETURN 

END 


CD 511) 
4 D 512) 
4 D 513) 
CD 514) 
CD 515) 
CD 516) 
CD 517) 
CD 518) 
CD 519) 
CD 520) 
CD 521) 
CD 522) 
CD 523) 
CD 524) 
CD 525) 
CD 526) 
CD 527) 
CO 528) 
ID 529) 
CD 530) 
CD 531) 
CD 532) 
CD 533) 
(D 534) 
CD 535) 
CD 536) 
CO 537) 
CO 538) 
CD 539) 
(D 540) 
CD 541) 
40 542) 
CD 543) 
CD 544) 
CD 545) 
CO 546) 
CD 547) 
CD 548) 
CD 549) 
CD 550) 
CD 551) 
CD 552) 
CD 553) 
CD 554) 
CD 555) 
CD 556) 
CD 557) 
CD 558) 
CD 559) 
CD 560) 
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APPENDIX D - Continued 

SURRClJT INF PH I 1 .1 

COMMON R SURO • AOVFRH,CHV< 2 ) • S 1 N XM 2 ) * CO SXH ( 2 ) • 

1 S I N XV ( 2 ) «Ct)SX V ( 2 ) * I HV90 t I CH I H* I CH IV* 

2 COORD < 3 1 # PH I | 2 ) fltJ.KfAHPRIM 

rsuro=sort< coord i i ) ** 2 *coord( 2 )**?+coordc 31**21 

AOVFRH=RSUBO*CHVI I HV90)-CO0RCC n*S INXHC ICHl H)*SINXVC ICHIV )+ 

1 AHPRIM*CnnRn(2)*SINXV( I CH I V ) *CGSXHC ICH I H ) + 

2 COORDC 3)*SINXH< I CH IH) *CC SXV C ICHIVI 
CAI I PARTI S(l«I«0'PTLin 

tf< i.fq.jigg rn ion 

r.AlL PARTLSI1 * J*0.PTL1 J) 

CAI L PARTL S C 3 • I *J*PTL2IJ) 
on TO 200 
100 PTL l J=PTL1 I 

CAI I P ARTL S < 2 • I.0.PTL2IJ) 

200 PHI! I ) = PTL2 I J/ AOVE RH— PTL1 1 *PTL 1 J/ AOVERH**? 

RFTURN 

FND 

SURROUT I NF PHI I1JK 

COMMON RSLJBO. A0VERH,CHV(2) • S I NXH I 2 ) • COSXH C 2 ) • 

1 SINXVC 21 .C0SXVC2) , IHV90* ICHIH. ICHI V, 

2 COOROI 3J«PHII2)fItJtK»AHPRIM 
IFIK.FO.UGO TO 100 

RSURO=SQRT( COOROI l )**2+C00RDC2 ) **2+ COORDC 3 1**21 

AOVE RH= RSIJB0*C HV ( I HV90) -COORD! l )*S INXHI ICHI H) *SINXV( ICHIVI* 

1 AHPR I M*COORD< 2)*SINXV( I CH I V > *COSXHI ICHIH)* 

2 COORD ( 3l*SlNXH( I CH I H ) *CC SX V ( ICHIVI 
100 CAM PARTLS(1*I*0*PTL1I1 

I F ( I.FO.JIGO TO 200 
CAI l PARTL SI1*J*0*PTL1J) 

GO TO 300 
200 PTL 1 J=PTl 1 I 

CAI I PARTLSI2* I • 0* P TL2 I J ) 

TFC I.FQ.KIGO TO 600 
GO TO 700 

300 IFM.FO.KIGO TO 500 
T F I J • FO . K ) GO TO 400 
CALL PARTLSII»K«0« PTL1K) 

CAI I PARTLSI3*J*K*PTL2JK) 

CALI PARTLSI3* I*K*PTL2IK) 

CAI I PARTI S 1 3 • I • J* PTL 2 I J ) 

CAM PARTLS I 6 * 0* 0« PTL I JK ) 

GO TO 800 
400 PTL 1 K^PTl 1 J 

CAI 1 PARTLS ( 2, J.0.PTL2JK) 

CAM PARTLSC3* I.K.PTL2IK) 

PTI 21 J=PTI 2 l K 

CALL PARTLSI5, J.I.PTLI JK) 

GO TO 800 
500 PTI 1 K^PTI 1 I 

CAI I PARTI S I 2 • MO* PTL2 IK I 
CAL l PARTL SI 3* J.K.PTL2JK) 

PTL 2 I J= PTL2 JK 

CAI I PARTLSI5* I*JtPTLIJK) 

GO TO 800 
600 PTL 1 K=PTLl I 

PTI 2IK=PTl 21 J 
PTI 2 JK=PTL2 IJ 
CAL L PARTL SI 4* I *0, PTL I JK ) 

Gfl TO 800 

700 CALL PARTLS(1.K*0,PTL1K) 

CAI I PARTL SI3.J.K*PTL2JK) 

PTL 2 I K— PTL2JK 

CALL PARTLSC5.J.K.PTLI JK) 

800 PHM2I = PTI I JK /AOVFRH— ( PTL1 I*PTL2JK+PTL 1J*PTL2IK 

1 +PTI IK*PTL2EJ)/A0VERH**2* 

2 2.0*PTLXI*PTL1J*PTL1K/ACVERH**3 

RFTURN 

FND 


(D 

561) 

(D 

562) 

(0 

563) 

(D 

564) 

ID 

565) 

(0 

566) 

ID 

567) 

(0 

568) 

(D 

569) 

ID 

570) 

to 

571) 

ID 

572) 

(D 

573) 

<D 

574) 

ID 

575) 

ID 

576) 

ID 

577) 

(D 

578) 

(D 

579) 

ID 

580) 

ID 

581) 

ID 

582) 

(D 

583) 

<D 

584) 

to 

585) 

to 

586) 

(D 

587) 

(D 

588) 

(D 

589) 

(0 

590) 

ID 

591) 

to 

592) 

(D 

593) 

(D 

594) 

ID 

595) 

(D 

596) 

(0 

597) 

(D 

598) 

(D 

599) 

ID 

600) 

(D 

601) 

(D 

602) 

(D 

603) 

to 

604) 

ID 

605) 

(D 

606) 

(D 

607) 

(D 

608) 

(D 

609) 

ID 

610) 

(D 

611) 

(D 

612) 

(D 

613) 

(D 

614) 

to 

615) 

ID 

616) 

CD 

617) 

(D 

618) 

CD 

619) 

ID 

620) 

CD 

621) 

(0 

622) 

CD 

623) 

CD 

624) 

CD 

625) 

CD 

626) 

(0 

627) 

to 

628) 

(0 

629) 
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APPENDIX D — Continued 


SUBROUTINE PARTLS( I ORDER, M , N , DER l V > 

COMMON RSUBO. AOVERH.CHVI 21 ,S INXHI 2 ) . COSXHI 2 1 . 

1 SINXV(2) »COSX V ( 2 J , I HV90 , ICHIH. ICHIV, 

2 COORD(3I.PH (2) ,I»J*K»AHPRIM 

C,n TO ( 100 . 200, 300, 400, 500,600) , I ORDER 

FIRST 0R0FR PARTI ATS WITH RESPECT TO COORDINATE M 
100 CO Ton 10.120. 1301. M 

1 10 DFR I V= ( COORD! 1 l*CHV< IHV90) ) /RSU80-S INXHI ICHIH1* 

1 S I N X V I ICHIV) 

GO TO 700 

120 DFR I V= ( COORO ( 2.) *CH VI IHV90) 1 /RSU80+AHPR I M* 

1 S I NX V I I CHI V ) *C0SXH ( I CHI H 1 

CO TO 700 

130 DFR I V= ( COORD! 3 I *CH V ( IHV90) l/RSUBO+SI NXH( ICHIH)* 

1 C0SXV( ICHIV) 

CO TO 700 

SFCOND ORDER PARTIAL WITH RESPECT TO COORDINATE M 

200 DFR I V= ( CHVI IHV90I/RSUB0I* ( 1 . 0- ( COORC ( M I **2 ) / 

1 <RSUBO**2>> 

CO TO 700 

SFCOND ORDFR PARTIAL WITH RESPECT TO COORDINATES M AND N 

300 DFR I V=-COORDIMI*COORD(N)*CHVI IHV90I /RSUB0**3 
CO TO 700 

THIRD ORDER PARTIAL WITH RESPECT TC COORDINATE M 

400 DFR I V= ( 3.0*COOROI M 1 *CHVI IHV90 I /RSUB0**3 I * 

1 ICnORDIM)**2/RSUBO**2-l.O) 

CO TO 700 

THIRD ORDER PARTIAL WITH RESPECT TO COORDINATES M AND N 

<500 DFR I V = ( COORD ( N 1 *CH VI IHV90 ) /RSUB0**3) *( 3. 0* 

1 COORO t MI**2/RSUB0**2— 1*0) 

CO TO 700 

THIRD ORDER PARTIAL WITH RESPECT TO THREE VARIABLES 

600 DERIV=3.0*C0nRDm*C00RD( 2)*COCRO( 3 ) *CHVII H V90 > / 

1 RSUB0**<5 

700 RETURN 
FND 


CD 

6301 

(0 

631) 

(0 

632) 

(D 

633) 

CO 

634) 

CO 

635) 

(0 

636) 

CD 

637) 

CO 

638) 

CO 

639) 

(0 

640) 

(0 

641) 

(0 

642) 

(0 

643) 

(0 

644) 

(D 

645) 

<0 

646) 

CD 

647) 

(0 

648) 

(D 

649) 

CD 

650) 

CD 

651) 

(D 

652) 

CD 

653) 

CD 

654) 

CD 

655) 

CD 

656) 

CD 

657) 

CD 

658) 

CD 

659) 

CD 

660) 

CD 

661) 

CD 

662) 

CD 

663) 

(D 

664) 

CD 

665) 

CD 

666) 

CD 

667) 

CD 

668) 

CD 

669) 

CD 

670) 

CD 

671 ) 

CD 

672) 

CD 

673) 

CD 

674) 

CD 

675) 

CD 

676) 
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APPENDIX D - Concluded 


THE FOLLOWING DATA CARDS WERE USED TO OBTAIN THE COMPUTED RESULTS PRESENTED 
IN TABLES II A.MO Ills 


COLUMN NUMBER 

0000000001111 nil 1122222222223333333333444*4444445555555555666666666677777777778 

123456789012345O789C 123456789012 34567 890123 45678901234567 8901 2345678901234567890 
SANGLES CHI H<1) =60.#CHIV (1 )=30.$ 


CLOSED TUNNEL CASE 

VELOCITIES 

AND DERIVATIVES 


2 1.000 1 .000 

1.000 

0.375 

C. 250 

0. 125 

2 2.000 0.75C 

2.000 

0.750 

0.250 

0.500 

2 1.000 1.000 

1.000 

0. 000 

0.000 

0.000 

$ANGLES CHIH(1) = 30. 

,C H I V ( 1 ) = 60. $ 




CLOSED TUNNEL CASE 

VELOCITIES 

AND DERIVATIVES 


2 1.000 l.COO 

1.000 

0.375 

0.125 

0.250 

2 0.800 1.500 

C.50C 

0. 375 

C. 250 

0.12 5 

2 1. 000 1.000 
■a 

1.000 

0.000 

0.000 

0.000 

3 

SANGLES CHI Fill =45 • 

, CH IV ( 1 ) = 45 • $ 




CLOSED TUNNEL CASE 

VELOCITI ES 

AND DERIVATIVES 


2 1.000 1.000 
a 

l.OOC 

0. 000 

0.000 

0. 000 

3 

SANGLES CHIH(1J=6C. 

,CHIV<1)=30.$ 




GROUND EFFECT CASE 

VELOCIT I ES 

AND DERIVATIVES 


2 2.000 0 . 75C 

2.000 

0. 750 

0.250 

0.500 


MORE TYPICALLY, WHEN MULTIPLE SKEW ANGLES ARE REQUIRED, THE FIRST CARD OF 
EACH GROUP MIGHT ASSUME A FORM SIMILAR TO THE FOLLOWING CARD: 

$ ANGLES CHI H(i> =50., 60., 70., 80. ,90 . , CH I V (1) =45 . , 60 . , 75 . , 90 . t 


THE SECOND CARD CF EACH GROUP MIGHT ASSUME ANY OF THE FOUR FOLLOWING FORMS: 


CLOSED TUNNEL CASE 
CLOSED TUNNEL CASE 
CLOSED TUNNEL CASE 
GROUND EFFECT CASE 


VELOCITIES ONLY 

VELOCITIES AND LONGITUDINAL DERIVATIVES ONLY 
VELOCITIES AND DERIVATIVES 
VELOCITIES AND DERIVATIVES 


COLUMN NUMBER 

CCOGOOOOO 1111 111 1 11222 222222 23333 33333344444444445555 5555 55666666666677 77 7777778 
123 45 67 89 01 2 34i> £>7 8901 2345 67 89012 34567 89012345678901 234567 8901 2 34567890 123 45 67890 
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to 


APPENDIX E 


FLOW CHART FOR FORTRAN PROGRAM OF APPENDIX D 


Program 5TABIL 
Function: Main Program 




















Program STAB I L - Concluded. 





CD 

CO 


APPENDIX E - Continued 


















CO 


Subroutine GRNEFF 

Function: Calculate interference factors in ground effect 



APPENDIX E - Continued 
















Subroutine WALEFF 


Function: Calculate interference factors in a closed wind-tunnel 



APPENDIX E — Continued 

























YES 


CO 

05 


Subroutine WALEFF. - Concluded. 



APPENDIX E - Continued 





























Subroutine PH I UK 

Function: Calculate interference gradients in free air 
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TABLE I 


VALUES OF ^ , AND THEIR DERIVATIVES FOR THE CENTER OF LIFT IN GROUND EFFECT 


I 


| For conciseness, the following abbreviations are used in this t able and in appendix C: 
= sinx H , Sy=sinx v , c H = cosx H , Cy=cosx v , C HV = Z 1 ' cos2 *H cos2 *V J 


A A ' 
h or T 




s y a z 
9 h 9 h 


j3_ 

y3 

a i 


jL 

-3 


3^ Ql 
h h 


TT h 


B # B f 


4 °i 


Values for use in determining - 


tan x v 1 

* -tan x v , 7 — , 1 

v tan x H J 


S H C V 


-2 s h s v 


2 s v c h 


S H C V U 2 c 2 C 2\ 
^(CHV - S H s v j 


C HV 2 


C HV 


CrV^ 


R 3„ _ 2 
S H S V C V 


CrV 2 


s h 2 s v c h c v 2 

C HV 2 


^^(chv 2 -h 2 V) 

C HV 


C HV 4 


3s h 3s v 2c v 3 


s H 2s V c H c V 2 />, 2 2- 2^ 

4 \ C HV - 3s H s V ) 

c hv* 


<x v : 


■900 X V , l] 


C HV + s v 
S H C V 


s H 3c V 3 / r , 2 2„ 2\ 

— 4“V C HV - 3s H s V ) 
C HV 


3^(°H 2 - 2s v 2c h 2) 

C HV 


-aV(s H 2 - 2sv 2 OH 2 ) 

C HV 


^%4 v 2 -2 s„v) 


C HV‘ 


C HV 


3Sh 4 b v 2c h c v 3 

CHV 4 ’ 


•^( C «V + •») 
^(Chv - *v) 


C HV 


a* 


-II 


*'|o,0,-2] 

2 ( C HV - s H c v) 


-s H s v 

" S V C H 

"( C HV “ s H c v) 


C HV 

2 


C HV 




**H 


S H C V 


-2s V c H 


°H C V U 2 „ 2 C 2' 

7“ “ 2\ C „V S H S V , 


C HV Z 


Chv 2 


Chv 


s h 3 s v c v 2 

chv 2 


s h 2 s v c h c v 2 

C HV 2 


3 2 

H S V C V 2 „ 2„ 2 

;74" 'i c Hv ' s « Sv 


C HV 


3s h s v c v 3 

c H v 4 


c H v 4 


“( C HV 2 ' 3 s H 2 s V 2 ) 


S H C V ( r 2 2. 2\ 

— 4“\^HV _ JS H S V ) 


4 S ^( s h 2 -W) 

C HV 




S H S V C V 


( s v 2 ” 2 s h 2 c v 2 ) 


OhV^v!( Sv 2. 2Sh 2 Cv 2 } 


C HV' 


3 s h 4 s v 2 c h c v 3 

C HV 4 


I x v = 90° |~ tan Xy> 

C HV + ®v 
S H C V 

-^(Chv 

'^ (C HV 
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TABLE II 


NUMERICAL DEMONSTRATION OF SYMMETRY OF INTERFERENCE FACTORS 


(a) In wind tunnel 


INTERFERENCE FACTORS FDR STABILITY WORK IN A CLOSED TUNNEL 
AT A Pn I NT NEAR A VANISHINGLY SHALL MCDEL 


GAMMA = 2.CC0 7FTA = 7.000 ETA = .750 

X/H = .750 Y/H = .750 Z/H = .500 


CHItHI = 60.00 CHI ( V > = 30. CO 

IX. XI IX. Y> (X, Zl IY.XI (Y,Y) (Y,Z) (Z.X) 17, Y) (7.ZI 
OFI T A ( - I -1.3957 -.B978 -.1094 -.5995 -.6524 2.C500 -.4285 -.7758 -1.7166 
DEI TAI-.-.XI — . 37 7 1 .8491 1.4477 1.19H9 .8414 .4786 -.7179 .5457 1.1973 
OF I TAI-.-.Y) .8491 1.1433 .512? .8414 3.0425 -1. 7695 .5457 -.6569 1.5429 
DFlTAI-.-./l 1.4477 .5172 -.766? .4786 -1.7695 -4.2414 1.1923 1.5429 1.3747 


(b) [n ground effect 


INTERFERENCE FACTORS FOR STABILITY HOF K IN GRCUND EFFFCT 
AT A POINT NFAR A VANISHINGLY SHALL HODFL 


X/H =' .750 Y/H = .750 Z/H = .600 


CHI I HI = 60.00 CHI ( V ) = 30.00 

IX. XI IX. V) ( X , 7 1 IY.XI IY.YI IY,Z> 17. X) IZ.YI (7.7) 

OF I T A I - I -. 1647 -.7271 -.1595 -.7771 -.0732 .4127 -.0153 -.1898 -.3646 

DFI TAI-.-.XI -.1851 .0712 .2345 .0712 .1183 .2463 -.2538 .0167 .0662 

OF l TAt-.-.Yl .071 7 .1 183 .2463 . 1103 .4851 -.0771 .0167 -. 150? .3173 

fJFI TAI-.-.7I .7345 .7463 .0667 .2463 -.0771 -.5563 .0662 .31 73 .4039 
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TABLE [(I 


NUMERICAL DEMONSTRATION OF INTERCHANGE EQUIVALENCES FOR A VANISHINGLY SMAa MODEL 


(a) Arbitrary point in arbitrary tunnel 


interference factors for stabiiity work in a closed tunnel 

AT A POINT NEAR A VANISHINGLY SMALL MCDEL 


GAMMA a 2.000 ZFTA = 2.000 ETA = .750 

X/H a .750 Y/H = .750 1 / H ' .500 


CHICHI = GO. 00 CHI (VI * 30.00 

IX. XI (X.YI I X • 7 ) ( Y. X ) (Y.YI lY.ZI IZ.Xl C7.YI 17. Zl 
OFLTAI-.-l -I. 3457 -.*97 8 -.1094 -.5993 -.6524 2.0500 -.4285 -.7258 -1.2166 
DFL TA(-.— «X I -.3771 .8491 1.4477 1.1989 .8414 .4786 -.7179 .5457 1. 1973 
0FITAC-.-.Y1 .8491 1.1433 .5122 .8414 3.0425 -1.7695 .5457 -.6569 1.5429 
DFITAI-.-.71 1.4477 .5172 -.7662 .4786 -1.7695 -4.2414 1.1923 1.5479 1.3747 


INTERFERENCE FACTORS FOR STABILITY WORK IN A CLCSFD TUNNEL 
AT A POINT NEAR A VAN I SH t NGL Y SMALL MCOEL 


GAMMA * .500 7FTA = .800 ETA = 1.500 

X/H * .375 Y/H = .250 2/ H = .125 


CHI I HI = 30.00 CHIIVI » 60.00 

IX. XI IX. Y) IX, 21 IY.X1 (Y.YI (Y,Z1 17. X) IZ.Y1 (Z.Z1 
0EITA4-.-I -1.3957 -.1094 -.8978 -.4785 -1.2166 -.7258 -.5998 2.0500 -.6524 
OFITAI-.-.XI -.7541 2.R855 1.6981 -1.4358 2.3847 1.0915 2.397R .9571 1.6828 
OELTAf-.-.Yl 7.RR5S -1.5325 1.0244 2.3847 7.7495 3.0858 .9571 -8.4878 -3.5389 
OF l TAI — .71 1.6981 1.0244 2.2866 1.0915 3.0858 -1.3137 1.6878 -3.5389 6.0849 


(b) Arbitrary point with centrally located model in square tunnel 


INTERFERENCE FACTORS FOR STABILITY WORK IN A CLOSEO TUNNEL 
AT A POINT NFAR A VANISHING! Y SMALL MCDEL 


GAMMA * 1.000 ZFTA - 1.000 ETA = 1.000 

X/H * .375 Y/H = .125 Z /H = .253 


CHI I H I = 30.00 CHIIVI « 60.00 

IX. XI IX. Y) IX«ZI (Y.X1 I Y ? Y ) tY.Z) (Z.X1 IZ.Yl (Z.ZI 
DEI TAI-,— I -.6061 -.4876 -.5072 .1588 -.7637 -.3555 -.0791 .2414 -.6279 
DEI TAI'-,-. XI -.8988 .0958 -.1753 -.4217 -. 1837 -.1564 .1357 .3045 .0363 
OFITAI-.-.YI .0958 .4388 .7309 -.1837 .6796 .7169 .3045 -.7749 .3787 
OFI TAI-,— .71 -.1753 .7309 .4599 -.1564 .7169 -.2584 .0363 .3707 .6392 


INTERFERENCE FACTORS FOR STABILITY WORK IN A CLOSED TUNNEL 
AT A POINT NEAR A VANISHINGLY SMALL MCOEL. 


GAMMA = 1,000 ZFTA = 1.000 ETA = 1.000 

X/H * .375 Y/H = .750 Z /H = .125 


CHI1HI = 60.00 CHIIVI * 30.00 

IX. XI IX. Y| I X »7 1 IY.X) (Y.YI ( Y, Z I (Z.X) (Z,Y| IZ.ZI 
DEI T A < - 1 -.6061 -.5072 -.4076 -.0791 -.6279 .2414 . 1588 -. 3555 -.7637 
OFITAI-.-.XI 89R 8 -.1 753 .0558 . 1357 .0363 . 3045 -.4212 -. 1564 -. 1837 
OFI TAI-,— . Yl -.1753 .4599 .7309 .0363 .6392 .3787 -.1564 -.2504 .7169 
OFI TAI-.-. 71 .0558 .7309 .4368 .3045 .3787 -.7749 -.1837 .7169 .6796 
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TABLE III. - CONCLUDED. 


(c) At centrally located model In square tunnel 


INTERFERENCE FACTORS FOR STABILITY WORK IN A CLCSEO TUNNEL 
AT A POINT NEAR A VAN l SH I NGl Y SHALL HCOEL 


GAMMA - 1.000 ZFTA = 1.000 ETA » I. 000 

X/H = 0.000 Y/H a 0.000 1 / H * 0.000 


r.HiiHi = AO. on cmm - 30.00 

IX. XI IX, Y1 IX.ZI IY.XI I Y » Y 1 IY.ZJ IZ.X) (Z.YI IZ.ZI 
OF I T A I - I —.099*5 -.5752 -.718* -.0*61 -.7899 .0299 .*210 -.2350 -.8231 
DFITA1-.-.X1 -1.09TB -.*978 -.5092 -.233* -.5315 .3050 -.1301 -.318* -.9078 
OFITAI-.-.Yl -.*978 .2308 .873? -.5315 .603* .82*9 -.318* -.5837 .5570 
DEI T&I-.-.7 I -.509? .8732 .8670 .3050 .82*9 -.3700 -.9078 .5570 .7138 


INTERFERENCE FACTORS FOR STABILITY WORK IN A CLCSEO TUNNEL 
AT A POINT NEAR A VANISHINGLY SHALL MCDEL 


GAMMA = 1.000 ZETA = 1.000 ETA * 1.000 

X/H = 0.000 Y/H = 0.000 1/ H * 0.000 


CHI I HI = 30.00 CHKVI * 60.00 

IX. XI IX. YJ IX.ZI ( Y. X I (Y,Y> IY.Z1 IZ.XI IZ.YI 17.71 
DFITAt-.-l -.0995 -.718* -.5752 .*210 -.8231 -.2350 -.0*61 .0299 -.7899 
OFITAI-.-.XI -1.0978 -.5092 -.*978 -.1301 -.9078 -.318* -.233* .3050 -.5315 
0EITAI-.-.Y1 -.509? .8670 .8732 -.9078 .7133 .5570 .3050 -.3700 .82*9 
DELTA!-. -.71 -.*97B .8732 .2308 -.318* .5570 -.5837 -.5315 .82*9 .603* 


Id) At centrally located model in square tunnel with equal vertical and horizontal skew angles 


INTERFERENCE FACTORS FOR STABILITY WORK IN A CLOSFO TUNNFL 
AT A POINT NFAR A VANISHINGLY SMALL MODEL 


GAMMA a L .000 ZETA = 1.000 ETA = 1.000 

X/H =* 0.000 Y/H = 0.000 Z /H = 0.000 


CHI1HI = *5.00 CHIIVI =« *5.00 

‘X.xi IX. Yl I X * 7 1 IY.XI IY.YI IY.ZI IZ.XI 17. Yl IZ.ZI 
DELTAI-.-l -.1369 -.5555 -.5555 .228* -.5826 -.2016 .228* -.2016 -.5826 
DEITAI-.-.XI -.98*7 -.3277 -.3277 -.07*0 -.*351 -.2219 -.02*0 -.2219 -.*351 
OFI TAI-.-.YI -.3777 .*92* .752? -.*351 .2152 .5286 -.2219 -.1912 .5286 
DFITAI-.-./l -.3777 .7522 .*97* -.2719 .5286 -.1912 -.*351 .5286 .2152 
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TABLE IV 


NUMERICAL DEMONSTRATION OF FAILURE OF INTERCHANGE EQUIVALENCES FOR A WING OF FINITE SPAN 


(a) Arbitrary wake skew angles 


INTfcRFERENCE FACTORS FOR STABILITY TESTS IN A CLOSED TUNNEL 
AVERAGE INTERFERENCE OVER A SWEPT WING OF FINITE SPAN 


ET A= 1.300 
StGMA= .753 


UNI FORK LOADING 
Z ET4 = 1.030 

AL PHA= 0.000 
BE T 4= 0.300 


GAMHA= 1.300 
L AM BDA= 0.000 


CH11H) = 60.00 CHIIVI = 30.00 


IX, X) 

DELTA!-,-) -.1799 
DELTA!-, -,XI -1.2S30 
0EITA1— »— ,Y) -.2477 
OELTA!-,-.Z) -.4959 


(X, Y) 
— .40B1 
-.2477 
.38 76 
.6169 


(X,Z) 
-.7790 
-. 4959 
.6169 
.8724 


( Y,X) 
. C81 9 
-.1467 
-. 7041 
-.1537 


I Y, Y) 
-.7885 
-.7141 
.07 52 
1.3358 


(Y, Z> 
-.1506 
-.1537 
1.0358 
.0715 


I Z , X > 
.3100 
-.9346 
-.0209 
-.6377 


(Z ,Y> 
- .0566 
-.0 229 
-.7291 
.1853 


(Z.Z) 

-.7749 

-.6877 

.1853 

.7637 


INTERFERENCE FACTORS FOR STABILITY TESTS IN A CLOSEO TUNNEL 
AVERAGE INTERFERENCE OVER A SWEPT WING OF FINITE SPAN 


FT A= 1.V9 
SIG9A= .742 


UNIFORM LOADING 
ZET A= 1.000 

ALPHA= 0.900 
RFTA= 0.900 


GAMM A= 1.000 

L AMBOA= o.oor 


CHIIHI = 3 u ,u CMltV) = ft.'.' 1 " 


OELTA! 

DELTA!-, -,X) 
DELTA!-, -,YI 
OELTA! Z) 


( X,X) 
-.'321 
-1 .6113 
-.0757 
-1.1591 


1 *,Y> 
- .0419 
-.9757 
1.4936 
1 . 559 ^ 


( X , Z ) 
-• fl" , 24 
-1 .1591 
1.5590 

.1182 


( Y , X ) 
.6348 
.9671 
-1.8371 
-.9156 


I Y » Y ) 
-1.0274 
-1.8871 
.9415 
1.1250 


<Y,Z» 

-.3644 

-.9156 

1.1250 

-1.0086 


<Z,X) 

-.2292 

-.8274 

1.0977 

-.7607 


! Z.YI 
. 1852 
1.0977 
-.8390 
1.0499 


(Z.Z) 

-.9873 

-.7607 

1.0499 

1.6664 


(b) Equal wake skew angles 


INTERFFRENCE FACTORS FOR STABILITY TESTS IN A CLOSEO TUNNEL 
AVERAGE INTERFERENCE OVER A SWEPT WING OF FINITE SPAN 


ET A= 1.000 
S IGM A= .750 


UNIFORM LOAOING 
ZETA= 1.09C 
ALPHA= 0.000 
8ETA= C.CO-'i 


G AMM A= l.COO 

L AMB DA= O.OCO 


CHIIHI =• 45.90 CH I ( V ) = 45. 'M* 


< X , X I 

OFl T A! - 1 -.1487 

OFL TA!-,-,X) -1 .4063 

OELTA!-,-, Y) - • 306 I 

OELTA!-,-, Z) -.*7*2 


f X,YI 
-.6373 
-.3261 
.7082 
.0598 


IX, Zl 
-.7278 
-.5702 
.9598 
.6981 


! Y , X I 
.3088 
-.0628 
-.8662 
-. 5662 


! Y , Y I 
-.7173 
-. 866 ? 
.1448 
. 92 84 


!Y,Z) 

-.3057 

-.5662 

.9284 

-.0319 


<Z,X) 

.1394 

-.0490 

.2382 

-.4162 


IZ.Y) 

-.0192 

.2382 

-.6664 

.3368 


<Z,Z) 
-.6728 
-.4162 
.3 368 
.7154 
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Figure 2 .- Wake in free air. 
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Figure 5.- Wake and images near intersection of simple ground plane and wall. 






Figure 7.- Wake and image with single sidewall. 
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Figure 8.- Wake and images near intersection of single sidewall and ground plane. 
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Effective vertical skew angle, X v , deg 


(a) 6 X( x and its gradients. 

Figure 13.- Interference factors as a function of x v for a vanishingly small model in ground effect and centered in a closed 
rectangular tunnel having a width-height ratio of 15 ,. x H = 90°. 
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Effective vertical skew angle ,X V , deg 

(g) 6 ZX and its gradients. 

Figure 13.- Continued. 
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0 10 20 30 40 50 60 70 80 90 

Effective vertical skew angle.Xy.deg 

(h) 6 Z y and its gradients. 

Figure 13.- Continued. 
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Effective vertical skew angle, X v , deg 



(a) 6 X x and its gradients. 

Figure 14.- Interference factors as a function of Xy for a vanishingly small model in ground effect and centered in a closed 
rectangular tunnel having a width-height ratio of 1.5. x H = 60°. 
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0 2 ' 


E f f ect i ve ver ti 


(b) 6x >y 
Figure 









I I I I I I 


Tunnel , y- 1.5 
Ground effect 






(e) 6y f y and its gradients. 
Figure 14.- Continued. 


Tunnel, y -\. 5 — 



3 30 40 50 60 70 


Effective vertical skew angle ,x V) deg 


(f) 6y z and its gradients. 
Figure 14.- Continued. 






20 30 40 50 60 70 80 9 < 

Effective vertical skew angle,X v ,deg 

(h) 6 z>y and its gradients. 


Figure 14.- Continued. 






Effective vertical skew angle, X v , deg 

(a) 6 x> x and its graaienis. 

Figure 15.- Interference factors as a function of x v tor a vanishingly small model in ground effect and centered in a closed 
rectangular tunnel having a width-height ratio of 1.5. x u = 30°. 









Effective vertical skew angle ,X v ,deg 

(c) 6 X z and its gradients. 

Figure 15.- Continued. 
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0 1 0 20 30 40 50 60 70 80 90 

Effective vertical skew a ng le , Xv,deg 


(d) 6y (X and its gradients. 
Figure 15.- Continued. 
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Effective vertical skew angle,X v ,deg 


(e) 6y y and its gradients. 


Figure 15.- Continued. 




Effective vertical skew angle, Xy,deg 

(f) 6y >z and its gradients. 

Figure 15.- Continued. 
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x,x,x' 


T unnel , y = 1.5 
Ground effect 


O I 1 1 1 I I ! 1 1 1 I 1 I 1 1 L i l I 

0 10 20 30 40 50 60 70 80 90 


Effective horizontal skew angle, X h, deg 

(a) 6 X x and its gradients. 

Figure 16.- Interference factors as a function of x H for a vanishingly small model in ground effect and centered in a closed 
rectangular tunnel having a width-height ratio of 1.5. Xy = 90°. 
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Figure 16.- Continued. 




Effective horizontal skew angle, X H ,deg 


(e) by y and its gradients. 
Figure 16.- Continued. 









I 



10 20 30 40 50 60 70 


Effective horizontal skew angle ,Xh, deg 

(i) 6 Z z and its gradients. 


Figure 16.- Concluded. 






s 



Effective horizontal skew angle, X h, deg 

(a) 6 X x and its gradients. 

Figure 17.- Interference factors as a function of x u f° r a vanishingly small model in ground effect and centered in a closed 

H 

rectangular tunnel having a width-height ratio of 1.5. Xy = 60°. 
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Tunnel , y- 1.5 
Ground effect 


1 



Effective horizontal skew angle .Xn.deg 

(b) S x> y and its gradients. 


Figure 17.- Continued. 






T unnel , y - 1.5 
Ground effect 



33 30 40 50 60 70 80 90 

Effective horizontal skew angle .Xn.deg 


(e) 6y y and its gradients. 
Figure 17.- Continued. 
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20 30 40 50 60 70 80 90 

Effective horizontal skew angle ,X h, deg 

(f) 6 y z and its gradients. 


— Tunnel , y = 1.5 
-- Ground effect 



20 30 40 50 60 70 80 

Effective horizontal skew angle, X h, deg 


(g) 6 Z x and its gradients. 
Figure 17.- Continued. 
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Effective horizontal skew angle, Xn, de 9 

(h) 6 Z y and its gradients. 

Figure 17.- Continued. 
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Effective horizontal skew angle,XH,deg 

(i) 6 ZjZ and its gradients. 

Figure 17.- Concluded. 
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Effective horizontal skew arigle.Xn.deg 

(a) 6 X x and its gradients. 

Figure 18.- Interference factors as a function of for a vanishingly small model in ground effect and centered in a closed 
rectangular tunnel having a width-height ratio of 1.5. x v = 30°. 




20 30 40 50 60 70 80 9C 

Effective horizontal skew angle, X H , deg 


(b) 6 Xf y and its gradients. 


Figure 18.- Continued. 



Effective horizontal skew angle, X h, deg 

(c) 6x 2 and its gradients. 

Figure 18.- Continued. 
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Effective horizontal skew angle, X h, deg 


(e) 6y f y and its gradients. 
Figure 18.- Continued. 
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Effective horizontal skew angle ,X h, deg 

(g) 6 ZX and its gradients. 

Figure 18.- Continued. 
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Effective horizontal skew angle,XH,deg 


(h) b Zf y and its gradients. 
Figure 18.- Continued. 
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Effective horizontal skew angle, X h, deg 

<i) 6 Z z and its gradients. 

Figure 18.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 19.- Distribution of interference factors along the longitudinal axis of the tunnel for a vanishingly small model centered 

in a closed rectangular tunnel having a width-height ratio of 1.5. x u = 90° ; x w = 90°. 
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(b) Caused by forces in the Y-direction. 
Figure 19.- Continued. 
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(c) Caused by forces in the Z -direction. 


Figure 19.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 20.- Distribution of interference factors along the longitudinal axis of the tunnel for a vanishingly small model centered 

in a closed rectangular tunnel having a width-height ratio of 1.5. x u = 90 °; X v , = 60°. 
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(a) Caused by forces in the X-directiori. 

Figure 21.- Distribution of interference factors along the longitudinal axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x. , = 90° • x*, = 30°. 
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(a) Caused by forces in the X-direction. 

Figure 22.- Distribution of interference factors along the longitudinal axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x H = 60°; xv = 90 °* 




(b) Caused by forces in the Y-direction. 
Figure 22.- Continued. 
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(c) Caused by forces in the Z -direction. 
Figure 22.- Concluded. 
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(a) Caused by forces in the X -direction. 

Figure 24.- Distribution of interference factors along the longitudinal axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x H = 60° ; x v = 30°. 
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(c) Caused by forces in the Z -direction. 
Figure 24.- Concluded. 




(a) Caused by forces in the X-direction. 

Figure 25.- Distribution of interference factors along the lateral axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x u = 90°; x = 90°. 
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(b) Caused by forces in the Y-direction. 
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(c) Caused by forces in the Z-direction. 
Figure 25.- Concluded. 
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Calculated 


From gradient 

Directly 
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(b) Caused by forces in the Y-direction. 


Figure 26.- Continued. 
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Calculated 

From gradient 
Directly 



(c) Caused by forces in the Z-direction. 
Figure 26.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 27.- Distribution of interference factors along the lateral axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x H = 90°; x v = 30°. 
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(c) Caused by forces in the Z-direction. 
Figure 27.- Concluded. 









(c) Caused by 1 
Figure 


Calculated 

From gradient 
Directly 




(a) Caused by forces in 


Figure 29.- Distribution of interference factors along the lateral axi 
in a closed rectangular tunnel having a width-heic 















(a) Caused by forces in the X-direction. 

Figure 30.- Distribution of interference factors over the lateral axis of the tunnel for a vanishingly small mode! centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. X H = 60°; x v = 30°. 
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(c) Caused by forces in the Z -direction. 
Figure 30.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 31.- Distribution of interference factors over the vertical axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x, . = 90°; Xw = 90°. 




(b) Caused by forces in the Y -direction. 


Figure 31.- Continued. 
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(c) Caused by forces in the Z-direction. 
Figure 31.- Concluded. 
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(c) Caused by forces in the Z -direction. 
Figure 32.- Concluded. 




(a) Caused by forces in the X-direction. 

Figure 33.- Distribution of interference factors over the vertical axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. X H = 90°; x v = 30°. 
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(a) Caused by forces in the X -direction. 

Figure 34.- Distribution of interference factors over the vertical axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. Xu = 60°; Xy = 90°. 




(c) Caused by forces in the Z -direction. 
Figure 34.- Concluded. 


(a) Caused by forces in the X-direction. 


Figure 35.- Distribution of interference factors over the vertical axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. = 60°; x v = 60°. 
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(b) Caused by forces in the Y-direction. 
Figure 35.- Continued. 




(a) Caused by forces in the X-direction. 

Figure 36.- Distribution of interference factors over the vertical axis of the tunnel for a vanishingly small model centered 
in a closed rectangular tunnel having a width-height ratio of 1.5. x H = 60° ; Xy " 30°. 



(b) Caused by forces in the Y-direction 
Figure 36.- Continued. 
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Figure 38.- Effect of span-width ratio o on gradients of interference factors for unswept wings 
centered in a closed tunnel having a width-height ratio y of 1.5. Xh = 90°. Note that 
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Figure 38.“ Continued. 
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Figure 38.- Continued. 
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Figure 39.- Continued. 
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Figure 40.- Continued. 
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Figure 40.- Continued. 
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Figure 40.- Continued. 
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Figure 40. ~ Continued. 







(a) Caused by forces in the X-direction. 

Figure 41.- Distribution of interference factors over the longitudinal axis of the tunnel for a uniformly loaded unswept wing centrally 
located and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 90°; Xv = 90 °- 
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(c) Caused by forces in the Z -direction. 
Figure 41.- Concluded. 
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(a) Caused by forces in the X -direction. 

Figure 42.- Distribution of interference factors over the longitudinal axis of the tunnel for a uniformly loaded unswept wing centrally 
located and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 90°; Xy = 60°. 
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(a) Caused by forces in the X-direction. 

Figure 43.- Distribution of interference factors over the longitudinal axis of the tunnel for a uniformly loaded unswept wing centrally 
located and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 90 °; Xy = 30°. 
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(c) Caused by forces in the Z-direction. 
Figure 43.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 44.- Distribution of interference factors over the longitudinal axis of the tunnel for a uniformly loaded unswept wing centrally 
located and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 60°; Xy = 90°. 
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(b) Caused by forces in the Y -direction. 
Figure 44.- Continued. 



(c) Caused by forces in the Z -direction. 
Figure 44.- Concluded. 
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(a) Caused by forces In the X-direction. 


Figure 45. Distribution of interference factors over the longitudinal axis of the tunnel for a uniformly loaded unswept wing centrally 
located and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 60°; Xy = 60°. 
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(b) Caused by forces in the Y-direction. 


Figure 45.- Continued, 
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(c) Caused by forces in the Z-direction. 
Figure 45.- Concluded. 


(a) Caused by forces in the X-direction. 


Figure 46.- Distribution of interference factors over the longitudinal axis of the tunnel for a uniformly loaded unswept wing centrally 
located and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 60°; Xv = 30°. 







(c) Caused by forces in the Z-direction, 
Figure 46.- Concluded. 


Calculated 

From centerline gradient 
From average gradient 
Directly 




Ca Iculated 

From centerline gradient 
From average gradient 

Directly 





I 

,y= 

— 

— 

— 


— 







- OX 

\ 

L 







V-s, 

<,z 



! 

V-8 

K,X 
















i 





i 














, j J i i i i i i i i i i i i 

-12 - 1.0 -.8 -.6 -4 -2 0 



(a) Caused by forces in the X -direction. 

Figure 47.- Spanwise distribution of interference factors for a uniformly loaded unswept wing centrally located and spanning 
half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 90°; Xy = 90°. 
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(b) Caused by forces in the Y-direction. 
Figure 47.- Continued. 
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Ca Iculated 
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(c) Caused by forces in the Z -direction. 
Figure 47.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 48.- Spanwise distribution of interference factors for a uniformly loaded unswept wing centrally located and spanning 
half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 90°; Xy = 60° - 





Ca Iculated 


From centerline gradient 
From average gradient 

Directly 




Ca Iculated 

From centerline gradient 
From average gradient 

Directly 



(a) Caused by forces in the X-direction. 

Figure 49.- Spanwise distribution of interference factors for a uniformly loaded unswept wing centrally located and spanning 
’ half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 90 °; Xv = 30 °- 


nil II I III Mil Hill I 


i mil 


I 1 1 Mil II 









288 


8 


(c) Caused by forces in the Z-direction. 



Figure 50.- Concluded. 



(a) Caused by forces in the X-direction. 


Figure 51. 


Spanwise distribution of interference factors for a uniformly loaded unswept wing centrally located and 
half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 60°; Xy = 60° ■ 
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(b) Caused by forces in the Y-direction 
Figure 51.- Continued. 
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(a) Caused by forces in the X-direction. 

Figure 52.- Spanwise distribution of interference factors for a uniformly loaded unswept wing centrally located and spanning 
half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 60°; Xy = 30°. 







(c) Caused by forces in the Z-direction. 
Figure 52.- Concluded. 
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(c) Caused by forces in the Z-direction. 
Figure 53.- Concluded. 
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<a) Caused by forces in the X -direction. 

Figure 54.- Distribution of interference factors over the vertical axis of the tunnel for a uniformly loaded unswept wing centrally Ixated 
and spanning half the width of a closed rectangular tunnel having a width-height ratio of 13. xh = 90 °; Xv = 60 °* 
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Figure 55.- Distribution of interference factors over the vertical axis of the tunnel for a uniformly loaded unswept wing centrally located 
and spanning half the width of a closed rectangular tunnel having a width-height ratio of 15. Xh “ 90 °; Xy * 30°. 





(b) Caused by forces in the Y-direction. 
Figure 55.- Continued. 
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(c) Caused by forces in the Z-direction. 
Figure 55.- Concluded. 
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(a) Caused by forces in the X-direction. 

Figure 56.- Distribution of interference factors over the vertical axis of the tunnel for a uniformly loaded unswept wing centrally located 
and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 6(P; Xy = 90°. 
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(a) Caused by forces in the X-direction. 

Figure 58.- Distribution of interference factors over the vertical axis of the tunnel for a uniformly loaded unswept wing centrally located 
and spanning half the width of a closed rectangular tunnel having a width-height ratio of 1.5. Xh = 60°; Xy = 30°. 
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(c) Caused by forces in the Z-direction. 
Figure 58.- Concluded. 
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Figure 61.- Concluded. 
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Figure 62.- Concluded. 
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